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FOREWORD 

This  report  is  part  of  a  study  whose  ultimate  objective  is  the  develop¬ 
ment  of  adequate  procedures  for  estimating  the  backsc altering  from  a  class 
of  rough  surfaces  at  oblique  angles  of  incidence.  In  order  to  avoid  a  com¬ 
plete  reliance  on  the  physical  optics  approximation  which  is  inaccurate  at 
angles  well  away  from  specular  and  does  not  provide  a  reliable  estimate  of  the 
de-polarizing  ***.  jt  of  the  surface,  the  approach  that  is  adopted  is  to  postulate 
a  simple  and  deterministic  base  surface  which  IS  itself  'rough*.  By  taking 
this  base  to  be  perfectly  conducting,  periodic  and  of  infinite  extent,  the  true 
surface  field  can  be  determined,  as  can  the  scattering  at  all  angles  of  in¬ 
cidence.  More  general  surface  configurations  such  as  might  be  appropriate 
to  the  sea  are  then  treated  by  a  modification  or  perturbation  of  the  field  on 
the  periodic  base. 

The  first  step  in  this  process  is  t"  investigate  the  scattering  of  a  plane 
wave  by  a  two-dimensional,  perfectly  conducting  periodic  surface.  This  is 
the  only  topic  considered  in  the  present  report. 
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ABSTRACT 


Numerical  procedures  are  developed  for  the  digital  solution  of  the 
integral  equations  for  the  currents  induced  in  a  perfectly  conducting, 
two  dimensional  periodic  surface  when  a  plane  electromagnetic  wave  is 
incident.  Data  are  obtained  for  both  the  surface  and  far  fields  for  a  variety 
of  sinusoidal  surfaces  for  plane  waves  of  either  polarization  at  oblique  as  well 
as  normal  incidence,  and  the  results  are  compared  with  the  predictions  of 
physical  optics. 
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I 

INTRODUCTION 

The  problem  of  the  scattering  of  an  electromagnetic  wave  by  a  periodic 
surface  is  one  which  is  easy  to  formulate  but  difficult  to  solve  rind  it  is  only 
recently,  with  the  advent  of  the  high  speed  computer,  that  reliable  results 
have  been  obtained. 

Prior  to  this  time,  a  variety  of  approximate  analytical  treatments  had 
appeared  in  the  literature,  most  of  them  being  based  on  the  approach  originated 
by  Lord  Rayleigh  (1878).  Assuming  the  surface  to  be  infinite  in  extent  and 
making  use  of  the  periodicity  in  (say)  the  coordinate  x,  Rayleigh  expanded  the 
scattered  field  in  discrete  spectrum  of  outgoing  plane  waves  ,  which  represen¬ 
tation  was  assumed  to  hold  right  down  to  the  surface.  Application  of  the 
boundary  condition  leads  to  a  single  equation,  valid  for  all  x,  from  which  to 
determine  the  complex  amplitudes  of  the  scattered  waves,  infinite  in  number. 

It  is  at  this  stage  that  approximations  must  be  made,  and  the  literature  is 
replete  with  methods  all  of  which  are  similar  in  character.  Thus,  for  a  per¬ 
fectly  conducting  sinusoidal  surface,  Rayleigh  (1878)  obtained  a  solution  by 
successive  approximation  based  on  the  initial  neglect  of  all  attenuated  waves; 
for  the  same  surface,  Tai  (  1948)  proceeded  via  an  orthogonal  mode  expansion, 
followed  by  matrix  truncation,  and  others  have  pursued  essentially  the  same 
path.  An  analogous  treatment  for  a  periodic  interface  between  two  homogeneous 
media  was  developed  by  Rayleigh  (1907),  and  applied  to  a  sinusoidal  profile  by 
Pavageau  (1963)  and  to  a  triangular  profile  by  Bousquet  (1963)  . 

At  best,  all  such  solutions  are  valid  only  for  corrugations  whose  height 
is  much  smaller  than  the  free  space  wavelength,  and  in  an  attempt  to  overcome 
this  restriction,  Meecham  (1956)  used  a  variational  method  to  find  the  angular 
distribution  of  scattered  energy  for  a  perfectly  conducting  grating.  The  scat¬ 
tered  field  was  represented  as  a  linear  combination  of  known  solutions  of  the 


1 


1363-6-T 


wave  equation  whose  coefficients  were  obtained  by  a  least  squares  fit  to  the 
boundary  condition  at  the  surface,  and  the  procedure  was  then  applied  to  a 
triangular  (or  saw  tooth)  profile.  A  somewhat  different  approach  was  taken 
by  Eckart  (1953)  and  Brekhovskikh  (1952)  who  both  used  variations  on  tht 
Kirchhoff  approximation,  and  by  Senior  (1959)  who  employed  physical  optics. 

In  either  case  the  determination  of  the  scattered  field  is  reduced  to  quadratures, 
and  Senior  showed  that  for  a  plane  wave  at  normal  incidence  on  a  sinusoidal 
grating,  the  physical  optics  integral  can  be  evaluated  exactly  to  give  the  com¬ 
plex  amplitudes  of  the  scattered  waves.  As  indicated  in  Appendix  A  of  this 
Report,  the  same  is  true  for  both  polarizations ,  and  for  oblique  incidence 
as  well  as  normal.  It  should  be  emphasized,  however,  that  the  solution  is 
still  approximate  by  virtue  of  the  postulated  surface  field  distribution,  and  the 
failures  of  the  physical  optics  estimate  of  the  surface  fields  are  examined  in 
Chapter  HI  of  this  Report. 

A  method  which  is  quite  distinct  from  all  of  the  above  was  developed  by 
Sivov  (1964)  who  used  conformal  transformation  and  a  consideration  of  the  Btatic 
limit  to  analyze  reflections  from  periodic  surfaces  with  shallow  and  deep  cor¬ 
rugations.  The  procedure  is  similar  to  that  recently  employed  by  Millar  (1969, 
1970)  to  investigate  the  inherent  limitation  of  Rayleigh's  method.  As  first  noted 
by  Lippmann  (1953),  it  is  not  in  general  valid  to  assume  that  the  expansion  of 
the  scattered  field  as  a  discrete  spectrum  of  outgoing  waves  alone  holds  over  the 
entire  scattering  surface,  and  this  fact  was  later  verified  by  Petit  and  Cadilhac 
(1964)  in  the  case  of  a  sinusoidal  grating,  hi  any  general  treatment  of  the  grating 
problem  it  is  therefore  necessary  to  allow  'ingoing'  waves  in  the  immediate 
vicinity  of  the  surface . 

Without  exception,  all  of  these  analytical  attempts  to  determine  the  fields 
scattered  by  periodic  surfaces  are  subject  to  approximation,  either  implicit  or 
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explicit,  and  it  is  only  with  the  use  of  high  speed  computers  that  have  permitted 
the  direct  digital  solution  of  integral  equations  that  reliable  results  have  been 
obtained. 

Most  of  the  initial  work  in  this  area  was  carried  out  by  the  French  inves¬ 
tigators  ,  such  as  Petit,  Cadilhac  and  Wirgin,  and  was  motivated  by  the  desire 
for  more  efficient  optical  diffraction.  In  his  early  papers.  Petit  (see,  for 
example,  1963)  followed  Rayleigh's  approach  in  expressing  the  scattered  field 
as  a  discrete  spectrum  of  outgoing  waves  alone,  leading  to  a  matrix  equation 
for  the  determination  of  the  spectral  amplitudes.  Since  the  matrix  was  then 
truncated  and  inverted  numerically,  it  will  be  appreciated  that  the  method  is 
no  more  than  a  digitization  of  that  originated  by  Rayleigh.  However,  in  later 
papers  (Petit,  1967),  the  Rayleig*  ssumption  was  circumvented  by  using  an 
integral  equation  formulation.  Series  expansions  were  adopted  for  the  incident 
and  scattered  fields  and  the  integral  equation  converted  to  a  matrix  equation 
which  was  solved  numerically.  Specific  results  were  obtained  for  plane  wave 
incidence  on  gratings  with  triangular  profiles,  and  the  efficiencies  computed. 

A  rather  different  approach  was  taken  by  Pavageau  (1967)  who  derived  the  in¬ 
tegral  equation  directly  in  terms  of  the  unknown  surface  current.  The  equation 
was  cast  in  the  form  of  a  nonhomogeneous  Fredholm  equation  of  the  second  kind 
and  solved  by  iteration. 

Methods  which  are  very  similar  to  that  which  we  shall  use  have  been  em¬ 
ployed  by  Neureuther  and  Zoiti  (1969)  and  by  Green  (1970).  The  former  con¬ 
sidered  scattering  by  periodic  structures,  either  dielectric  or  perfectly  con¬ 
ducting.  The  integral  equation  was  obtained  from  Green's  theorem  and  the 
modified  Green's  function  expressed  in  either  of  two  ways  depending  on  the 
parameters  of  the  surface.  The  first  (space  harmonic)  representation  is 
analogous  to  that  used  by  Petit;  the  second  consists  of  an  infinite  series  of 
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Hanks 1  functions,  &ud  was  computed  using  a  Mellin  transform  and  an  asymptotic 
comparison  scheme.  Numerical  data  were  obtained  for  sinusoidal  profiles. 
Green  0970)  also  used  the  spacs  harmonic  representation,  but  improved  itc 
convergence  by  summation  techniques.  Data  were  presented  for  the  surface 
fields  and  diffraction  efficiencies  of  perfectly  conducting  gratings  with  tri¬ 
angular  profiles. 

The  present  work  also  employs  the  numerical  approach,  and  is  concerned 
with  the  scattering  of  electromagnetic  waves  by  infinite,  perfectly  conducting, 
two  dimensional  periodic  surfaces  of  arbitrary  but  continuous  profile.  Plane 
wave  incidence  is  assumed,  with  either  £  or  H  polarization,  and  both  normal 
and  oblique  incidence  are  considered.  A  representation  of  the  Green':*  function 
is  employed  that  is  similar  to  that  used  by  Green  0970),  and  the  convergence 
improved  still  further  by  subtracting  the  de  terms.  This  has  the  added  advantage 
of  making  explicit  the  behavior  of  the  Green's  function  in  a  neighborhood  of 
its  singularity.  Digital  programs  have  been  written  for  computing  the  surface 
and  scattered  fields,  and  the  program  listings  are  included  in  Appendix  B. 

Data  are  presented  for  a  variety  of  sinusoidal  profiles,  and  the  results  compared 
with  the  physical  optics  predictions . 
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n 

FORMULATION  OF  THE  INTEGRAL  EQUATION 
FOR  NUMERICAL  SOLUTION 


We  consider  here  an  infinite,  perfectly  conducting  periodic  surface  illuminated 
by  a  plane  electromagnetic  wave.  Since  the  surface  is  assumed  two  dimensional 
in  the  sense  of  being  independent  of  a  Cartesian  coordinate  z,  the  entire  problem 
is  two  dimensional,  and  the  most  general  solution  can  be  deduced  from  the  par¬ 
ticular  solutions  appropriate  to  incident  plane  waves  having  either  E  or  H  in  the 
z  direction,  i.e.  parallel  to  the  corrugations.  In  either  case,  the  problem  is 
essentially  scalar. 

2. 1  Formulation 

It  is  convenient  to  develop  first  the  integral  equation  in  the  somewhat  simpler 
case  of  a  scattering  surface  of  finite  extent.  Let  S(see  Fig.  2-1)  be  this  surface, 
and  surround  it  by  another  closed  surface  SR  .  Let  0(r)  and  G(rjr')  be  two 
scalar  functions  which  are  continuous,  together  with  their  first  and  second 
derivatives,  on  S  and  and  throughout  the  volume  V  enclosed  by  them.  Assume  , 
moreover,  that  ijj( r)  satisfies  the  homogeneous  wave  equation 

(V2-Hc^(r)=0  (2.1) 

inside  and  on  the  boundaries  of  V,  whereas 

(  V2  +  k2)  G(r|r’)  =  -  6  (r-r»)  .  (2. 2) 

Applying  Green’s  theorem  to  the  volume  V,  we  obtain 
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/ 


G(r|r')  ^(r‘)-^(r')  g|,  G(rjr*) 


dS* 


■ft 

V 

=  i//(r) 


G  (r  |r»)  (V24k2)(//(r ')-0(r')  ( V2-*2)G  (r  |  r ') 


dV' 


(2.3) 


by  virtue  of  the  properties  of  the  delta  function.  We  now  identify  ip  (r)  with  a 
total  field 


^t0t  (l)  =  0*(l)  +  ^8(r) 


where  <//*(r)  is  an  incident  plane  wave  originated  by  a  source  at  infinity  (and 
therefore  outside  V),  and  08(r)  is  the  field  scattered  by  the  surface  S  .  Since 
ipB(r)  must  satisfy  a  radiation  condition  at  infinity,  its  contribution  to  the  in¬ 
tegral  over  3^  decreases  to  zero  as  the  surface  recedes  to  infinity,  whereas 
\pHr)  contributes  itself.  Equation  (2.3)  then  becomes 


.tot. 
ifj  (r)= 


[r)r  ^(r)  +  J  ^G(r|rr 


3  n 


9  ‘,  ^(r1)  -^(r1)^— ,  G(r|r') 


dS' 


(2.4) 


and  in  spite  of  the  assumptions  of  an  incident  plane  wave  and  a  surface  of  finite 
extent,  Eq.  (2.4)  is  also  valid  for  an  arbitrary  Incident  field  and  for  a  surface  S 
extending  to  infinity.  In  the  latter  case,  however,  the  proof  is  by  no  means 
trivial  (Jones,  1952). 

The  particular  situation  of  concern  to  us  is  that  in  which  S  is  doubly  infinite 
and  divides  apace  into  two  regions.  It  is  then  sufficient  to  integrate  over  only 
the  upper  ('illuminated'  )  side  of  the  surface.  We  also  assume  that  S  is  independent 
of  the  coordinate  z,  and  this  allows  us  to  distinguish  two  particular  cases  according 
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to  the  polarization  of  the  incident  field.  Mks  units  are  employed  and  a  time 
factor  e^4  suppressed  throughout. 

2.1.1  E  (or  horizontal)  Polarization 

If  E*  =  E*  z  .  the  electric  vector  in  the  scattered  field  will  also  be  confined 
-  z 

to  the  z  direction,  and  we  can  make  the  Identification 


<Hr)  =  Ez(r)  .  (2.5) 

Since  the  normal  derivative  of  Etot  is  related  to  the  surface  current  density  K 

z  z 

by  the  equation 


!;•  <r  & ")UM  KI<£,)  • 


(2.6) 


where  Mis  the  permeability  of  free  space,  Eq.  (2.4)  can  be  written  as 

e‘°V>  =  E‘z <£> / (G'l|l'»Kz(r'»+  j~E^ot(r')  A,  G(r|  r')  j  OS'  . 

s  (2.7) 


At  the  perfectly  conducting  surface  S  the  boundary  condition  is 

EtOt(r)=0  (2.8) 

z  — 

and  hence,  on  allowing  r  to  approach  the  surface,  Eq.  (2.7)  gives  rise  to  the 
integral  equation 
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J°fe 


Gfelr'W^W  '  jijE*  © 


(2.9) 


For  a  two-dimensional  problem  the  free  space  Green's  function  is 


G(r|r')  =  -  i  H^(k  |r-r'|  ) 


(2. 10) 


where  '  is  the  zero  order  Hankel  function  of  the  second  kind .  The  final 
integral  equation  is  therefore 


J Kz(r,)H®(k|r-r'|  Wl*  =^£*0 


(2. 11) 


from  which  K  has  to  be  determined, 
z 

2.1.2  H  (or  vertical)  Polarization 


1  _  J  * 


If  H*  -  H*  z  ,  the  total  magnetic  field  is  likewise  in  the  z  direction,  and  on 

2 

making  the  Identification 


*<£)  =  H‘ot(r)  . 


(2. 12) 


Eq.(2.4)  becomes 


Hz  fe)+J  S’.CMrlr')  J  d8'  . 

S  1  (2. 13) 


The  boundary  condition  at  the  perfectly  conducting  surface  S  is 
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a_ 

3n' 


(r)  =  0 


(2.14) 


and  on  allowing  r  to  approach  the  surface,  Maue  (1949)  has  shown  that  (2. 13) 
reduces  to 


S 


Htz0t(r,)~G(r|r,)  dS*  . 


(2. 15) 


The  quantity  H  is  the  induced  current  density,  and  since  this  current  flows 
z 

tangential  to  the  surface  in  a  plane  perpendicular  to  the  z  direction,  we  can  write 


K, (r)  =  Ht0t(r)  t 
— t  ~  z  - 


(2. 16) 


where  t  is  a  unit  vector  tangential  to  S.  Inserting  finally  the  expression  (2. 10) 
for  the  free  space  Green's  function,  we  have 


JK  tte'sX* 


Qc  |  r-r'|  )  df '  =  4J 


h‘ (r)  -  i  Kt  <r> 


(2. 17) 


i.e. 


J Kt(r)H(J2)(k|r-r'|  )cos(n  ,|  r-r' j  )d<«  -  4j|hVj  -  ^(rlj  .  (2. 18) 
S 


This  Is  an  Integral  equation  from  which  to  determine  . 

2. 2  Reduction  of  the  Integral  Equations  fox  Periodic  Surfaces 

We  now  make  use  of  the  fact  that  the  ncident  field  ia  a  plane  wave  and  that 
the  surface  y  =  f(x)  ia  periodic  with  period  d,  i.e. 
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f(x  +  md)  =  f(x),  m  =+  1,  +  2,  +3,...  .  (2.19) 

As  before,  it  is  convenient  to  consider  separately  the  cases  of  E  and  H  polari¬ 
zation. 

2.2.1  E  Polarization 
Let  us  assume  that 

El  (r)=e‘ik(x!,Ine‘yco89)  (2.20) 

Z  “ 

where  0  is  the  angle  of  incidence  with  respect  to  the  normal  to  the  mean  surface. 
The  integral  equation  is  as  shown  in  Eq.  (2.11),  but  since  the  integration  extends 
from  -  oo  to  a>,  this  is  not  appropriate  for  a  numerical  solution.  However,  by 
invoking  the  periodic  property  of  the  surface,  the  integral  can  be  reduced  to  one 
over  a  single  period  alone  at  the  expense  of  a  mere  complicated  form  for  the 
kernel. 

From  Floquet's  theorem,  we  have  that 

K  (r  +  mdx)  =  K  (r)  e“jkmd  8ln9  ,  m  =  ±  1,±2,  +3,  . . .  (2.21) 

z  —  z  ~ 

Moreover, 

di*  =  M+ff(x*)J2  dx* 

where  f’(x)  is  the  derivative  of  f  with  respect  to  x.  This  allows  us  to  express 
the  integral  as  one  along  the  x'  axis,  and  since 

|  r-r ' |  =  , 
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we  obtain 


f  i tffr*51  rad)2+(y-y,>21  ^e"Jkmdain0  ^l+{f‘(x')}^kz(Jt')dx' 


4  -jk(xsin0-ycos0) 
uyt 


(2.22) 


valid  for  0  <  x  <  d  with  y  =  f(x),  y*  =  f(x'). 

In  arriving  at  (2.22)  we  have  implicitly  assumed  that  interchanging  the  order 
of  integration  and  summation  is  valid  (a  fact  which  is  by  no  means  obvious),  but 
even  so  Eq.  (2.21)  is  not  a  very  promising  equation  for  numerical  purposes 
because  of  the  extremely  poor  convergence  of  the  infinite  series  constituting  the 
kernel.  To  rectify  this  situation,  consider 


^  (k^x-x'-md^+fy-y'?  V 


-jkmdsind 


mr-oo 


(2.23) 


Using  the  Fourier  integral  representation  of  the  Hankel  function,  the  Poisson 
summation  formula  (see,  for  example,  Morse  and  Feahbach,  1953)  applied  to 
(2. 23)  gives 

-j(~2.+k  sin0)(x-x')  -j  |y-yH  X 

d  e  ",  (2.24) 


where 


X 

m 


—  . . . ~i 

+  k  sin#)2 


(2.25) 


and  the  chosen  branch  of  the  square  root  Is  that  for  which 


ImX  <  0 
m  - 
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in  order  to  satisfy  the  radiation  condition. 

The  expression  (2. 23)  for  P^has  been  used  by  Yen  (1962)  and  Green  (1970),  and 
is  equivalent  to  that  obtained  by  Neureuther  and  Zakl  (1969)  using  a  space  har¬ 
monic  expansion.  The  features  that  should  be  noted  are  that  if  yf  y'  the  suc¬ 
cessive  terms  decrease  exponentially  with  increasing  I  ml  for  all 


m  +  “  sine 
2  * 


and  that  in  addition  there  is  an  algebraic  decrease  (proportional  to  m"  *)  provided 

by  the  factor  X  in  the  denominator, 
m 

We  can  produce  still  a  further  improvement  in  the  convergence  properties 
of  the  series  by  separating  out  the  zero-frequency  (k=0)  terms,  and  this  has  the 
added  advantage  of  making  explicit  the  behavior  of  P^in  the  vicinity  of  the  singular 
point  x'=x,  y,=y  .  Putting  k  =  0  in  (2.24),  we  have 


|k=0 
and  since 
oo 


2 

d 


2  JSE 

ra=-oo  Jlralx 


1  -J^(x-x')  -^|m(y-y*)| 


2  it 


(2.25) 


2 

m=l 

■1 


2irJ 


log 


-j  tt  <*-**>-  t  I  y-y’1  ) 


1-e 


f  j~(x-x()  -  J-  |y-y'| 
1-e 


m=-oo 


(see,  for  example,  Collin,  1960),  provided  y1  f=y  ,  it  follows  that 
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i 


P,|  =  (m=0  term)-  J  logi 

lk®0  * 


-  ^  [ly-y'l  tJC'-x*)] 

t 

-  j  [ly-y'l  -)0t-x»)j  | 


(2.  27) 


Hence,  by  subtracting  the  zero  frequency  value  of  eaoh  term  of  (2. 24)  and  then 
adding  (2.27),  the  integral  equation  (2.21)  takes  the  form 


d 

I 


jkx’sinfi  / — - - ~~i 

G  (x.y;  x’,y')K  (x’)e  V1+  W))2  dx’  =  ~  ejkycose 

1  z  i  J  uu  * 


0<x<d  ,  (2.28) 


where 


G.<x.y j  x-.y')  .Ti_B-!k  ly-y'l 
1  kcosS 


kCO80 


l-e 


-—Oy-y’i  -]<*-*’>) 


*  1 


a,  -J ly-y'l ^ 


1 


2mjr 


m  d 

+  "5m*r 


ly-y’l 


m 


J* 


A  -Jiy-n*;  -t1  ly-y’l) 

Z*  d  h. - +  y 


+  4<e 

m-1 


2m* 


m 


J 

(2.29) 
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with 


We  observe  that  the  Green's  function  has  now  been  expressed  as  the  sum 
of  a  term  involving  explicitly  the  amplitude  and  the  angle  of  incidence,  a 
logarithmic  function  representing  the  true  behavior  of  G,  in  a  neighborhood 
of  its  singularity,  and  two  series  which  are  themselves  differences  of  two 
convergent  series. 

2.2.2  H  Polarization 

The  incident  field  in  this  case  is  taken  to  be 


H*  (r)  -  g-Jkfr  8in0-V  cosfi) 


(2.31) 


and  the  reduction  of  the  integral  equation  (2.  '7)  to  a  form  suitable  for 
numerical  solution  proceeds  in  much  the  same  way  as  for  E  polarization. 
Since 


Eq.  (2.  17)  can  be  written  as 

P2Kt(x')dx'*  4jf  H*  (r)-  ~  K^r)  1  (3. 32) 
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where 


/l+|f'(x*)}  2  —  ^  H*2) 


e 


-jkrodaiae 


m=-a> 


2  X  jf'Cx'K  IT  \  x 

ra=-CD  „ 

”j{”^+kslne)(x-xH  -j  jy-y'j 


(2.33) 


in  which  sgn(T)  is  the  signum  function  and  ia  given  in  (2.25).  The  zero  fre¬ 
quency  (k=0)  limit  of  Pg  is 

1  -  y  |m(y-y’)l 


2!k=0  d 


m=-oo 


2  1  agato"y'f 


d  d 

e  e 


2  (  )  2{  )  T  'j2*T 

=  ^|f,(x,H  sgn(y-y’)  1  f’fa'Hj sgafy-y’) j  2L  e 

m=l 


2 

+a 


-f*Cx*)+j  sgnfy-y’) 


^  Jx{^ ly-y's) 


m=l 


and  since 
oo 


V  imT  i  ,  i  *T 

L  e  '  ‘  2  +2  cot  2 


m=l 
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we  have 


k=0 


jj  f*(x*)+  j  |  f'CxO+j  sgn(y-y')|  cot  j-  ~  {x-x')+j|  |y-y'| 

|“f'Cx')+ j  sgn(y-y‘) | cot  (x-x'Rjrj  !y-y'|  |  .  (2.34) 


Hence,  by  subtracting  the  zero-frequency  value  of  each  term  of  (2. 33)  and  then 
adding  (2.34);  the  integral  equation  (2.32)  takes  the  form 


( 


G^.yjx'.y'Jl^tx'Jdx'  =  jd|2H^(x)-Kt&)  j  ,  0<x<d,  (2.35) 


where 


G2 fe,y?t,.y,)=  -j  egu(y-y'H  |f  (x*)tan9-sgn(y-y')}  Hkcosfl  ^y_y  ^ 

-3  ly-y'lx* 


®  Vx') 

_je-jksine(x-x«)  J  e  3  d  *  x  3 

ra=l 


X+ 

m 


f,(x')(^--H«sin0)-X+  sgn(y-y')l 
a  m  j 


2m7r 


iy-y’ 


-e 


Ijjf'fx'J-sgnty-y')  j 


.  -jksin0(x-x')  *  d 

-je  ?  e 


cd  42m£/v_x,)  -j|y-y*lx^ 


2 

m=l 


( 


2m7T 


2  m  7T, 
d 


f'(xO(--~j~+-ksin0)  -Xmsgn(y-y'^e  ^  jjf'fcO-tegnfy-y')}  J 


-|Jf(x,)+jsgn(y-y,)J  cot  ^  (x-x')+|  ^|y-y'l  j 

■i  { -f'Oc'M  sgn(y-y‘)j  cot  j  ^(x-x')+)  J  |  y-y*  |  J  .  (2. 36) 
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In  spite  of  tbs  obvious  complexity  of  the  Green's  function  Gg.  its  form  is 
directly  analogous  to  in  consisting  of  certain  explicit  terms,  a  pair  of  cotangent 
functions  representing  the  true  behavior  of  Ggin  a  neighborhood  of  its  singularity, 
and  two  series  which  are  themselves  differences  of  two  convergent  series.  Once 
again,  therefore,  the  singularity  of  the  Green's  function  has  been  separated  out. 

We  further  note  that  the  integral  equations  (2.28)  and  (2.35)  have  been 
derived  without  approximation  .  hi  consequence,  the  formulation  so  far  is  exact 
for  any  two  dimensional  perfectly  conducting  periodic  surface  which  is  smooth  in 
the  sense  of  having  a  continuous  first  derivative. 

2.3  Numerical  Procedures 

Methods  for  the  numerical  solution  of  integral  equations  have  been  extensively 
discussed  in  the  literature  (see,  for  example,  Harrington,  1968).  The  general 
procedure  consists  of  reducing  the  equation  to  a  finite  set  of  algebraic  equations, 
i.e.  to  a  matrix  equation,  and  can  be  illustrated  by  considering 

G(x,x')K(x')dx'  =  F(x)  ,  0  <  x  <  d  .  (2. 37) 


We  assume  that  tlie  unknown  function  K(x')  can  be  expanded  in  terms  of  linearly 
independent  base  functions  ^(x1)  such  that 
N 

K(x')£  S  a  0  (x')  (2.38) 

dart  Q  n 

n=l 


where  the  aQ  are  the  associated  constants.  Substitution  of  (2.38)  into  (2. 37) 
gives 


G(x,x')0  (x')dx'  =  F(x) 
n 


(2. 39) 
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and  the  solution  of  the  integral  equation  has  now  been  reduced  to  the  determination 
of  the  constants,  aQ,  n=l,2, . . .  N  . 

There  are  several  possible  ways  of  finding  the  an>  e.g.  least  squares  fit, 
Galerkin's  method,  and  the  collocation  method,  and  it  is  the  last  of  these  that 
we  shall  employ. 

The  collocation  method  converts  (2.39)  into  a  system  of  N  linear  equations 
by  forcing  the  two  sides  of  (2. 39)  to  be  equal  at  N  sampling  points  in  the  interval 
(0,  d).  This  is  simply  a  point  matching  procedure  and  results  in  the  matrix 


equation 


G(x  ,x')0  (x')dx'  =  F(x  )  , 
m  n  m 


0<x  <  d,  m=l,2, 

—  m  — 


.N  . 


(2.40) 


There  now  remains  the  problem  of  choosing  the  base  functions  0q(x*)  ,  and  here 

again  there  are  several  possible  choices,  e.g.  rectangular,  quadratic  and 

sinusoidal.  By  appropriate  choice,  we  can  economize  in  the  number  of  sampling 

points  required  for  an  'accurate'  approximation  to  the  solution  K(x),  and 

experience  has  shown  that  a  rectangular  function  is  not  in  general  a  good  choice, 

whereas  sinusoidal  interpolation  often  works  rather  well. 

The  particular  form  of  sinusoidal  interpolation  that  we  have  adopted  is 

predicated  on  the  use  of  sampling  points  which  are  uniformly  distributed  in 

0  <  x'  <  d  .  The  range  of  integration  is  therefore  broken  up  into  N  increments 

of  length  A  =  d/N  (see  Fig.  202).  Furthermore,  let  x'  be  the  midpoint  of  the 

n'th  cell,  i.e.  x'  =  (n-x)A,  and  let  Ax'  denote  the  interval  x'--*-  <  x'  <x'  +  4r  . 

nz  n  n  z  —  —  n  2 

We  assume  that 


K(x')  =  { 


A  +B  sin  k(x'-x  )+C  cosk(x'-x  ) 
n  n  n  n  n 


V.  0 


,  if  x'£  Ax' 
n 

otherwise 


(2.41) 
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FIG.  2-2:  flluetraUoa  of  the  inter poUtion  procedure. 
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the  periodicity  of  the  ^oblem  must  be  used  to  determine  the  constants  K  or 
K_TJ,  required  in  (2.  *3) 

The  ’bove  procedure  is  immediately  applicable  to  the  integral  equations 

for  E  and  H  polarizations  on  inserting  the  appropriate  values  for  the  Green's 

function  and  the  forcing  function,  but  a  brief  comment  is  desirable  concerning 

the  treatment  of  the  singular  cell  .  The  Green's  functions  of  concern  to  us  are 

singular  when  x  -  x'  ,  the  singularity  being  logarithmic  for  E  polarization,  and  a 

first  order  pole  for  H  polarization,  and  it  is  therefore  necessary  to  modify  the 

numerical  scheme  when  x  =  x  .  In  line  with  the  usual  practice,  we  divide  the 

m  n 

singular  cell  into  three  portions: 


Vi-Vi1' 


vi> 


with  0  <  c  <  A .  The  first  and  third  segments  are  handled  by  the  standard 
numerical  technique,  whereas  the  central  portion  is  treated  analytically  by  means 
of  a  limiting  process  (Andreasen,  1964). 
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m 

APPLICATION  OF  THE  NUMERICAL  TECHNIQUE 
AND  PRESENTATION  OF  DATA 


3. 1  Discussion 

Although  the  only  data  that  we  shall  present  are  for  sinusoidally  cor¬ 
rugated  surfaces,  it  is  well  to  remember  that  the  formulation  given  in  Chapter 
II  is  quite  general  and  appropriate  to  any  perfectly  conducting  .two-dimensional 
periodic  surface.  It  is  not  even  necessary  that  the  surface  have  continuous  first 
derivative,  and  providing  care  is  taken  of  any  singularities  at  sharp  edges,  the 
procedure  is  applicable  to  surfaces  whose  first  derivative  is  only  piecewise 
continuous  as  in  the  case  of  a  saw-tooth  profile.  The  particular  numerical 
scheme  that  has  been  adopted  is  especially  suitable  for  surfaces  having 
2jra/d»  1,  where  a  is  the  peak  deviation  of  the  surface  from  its  mean,  and 
d  is  the  period.  In  the  infinite  series  representation  of  the  Green’s  function 
(see  Eqs.  (2.29)  and  (2.  36) ),  the  first  few  terms  (for  which  1  sin  flj<  1  ) 
are  oscillatory,  but  for  larger  m  the  crms  decrease  exponentially  at  a  rate 
which  is  ultimately  determined  by  the  quantity  2?ra/d.  Whereas  the  earlier  terms 
correspond  to  propagating  modes  in  the  scattered  field,  the  later  opes  correspond 
to  evanescent  modes,  and  it  is  only  necessary  to  include  a  few  of  the  latter  in 
order  to  achieve  adequate  accuracy.  It  is  therefore  apparent  that  for  surfaces  of 
relatively  small  period  (  d£  A)  the  infinite  series  can  be  approximated  by  the 
first  few  terms  alone. 

The  unknown  function  in  the  integral  equation  is  the  current  induced  in  the 
surface,  and  this  is  the  quantity  that  is  computed  initially.  Knowing  the  surface 
field  it  is,  of  course,  a  trivial  matter  to  determine  the  scattered  field,  and  the 
means  of  obtaining  this  in  the  form  that  is  most  convenient  for  our  purposes  is 
described  in  Section  3.2. 
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As  noted  earlier,  all  of  the  numerical  results  given  in  this  report  are  for 
sinusoidal  surfaces  whose  profile  is  taken  to  be 

y  =  a  cos  .  (3. 1) 


Several  different  combinations  of  d/X,  a/X  and  incidence  angle  are  considered, 
and  the  distributions  of  surface  and  scattered  fields  computed  .  Comparisons 
with  the  predictions  of  physical  optics  are  also  given. 

3.2  The  Scattered  Field 

When  a  plane  wave  is  incident  on  a  periodic  surface,  the  scattered  field 
can  be  represented  as  an  angular  spectrum  of  plane  waves,  which  spectrum  is 
discrete  by  virtue  of  the  periodic  nature  of  the  boundary  condition  at  the  sur¬ 
face.  Each  of  the  infinity  of  waves  making  up  the  spectrum  has  associated  with 
it  a  diffraction  angle  which  may  be  real  or  complex  and  is  determined  by  the 
grating  law.  Whereas  the  amplitude  of  the  wave  is  a  function  of  the  profile  size 
and  shape,  and  the  directions  of  incidence  and  diffraction,  the  diffraction  angle 
depends  only  on  the  value  of  d/X  and  the  direction  of  incidence. 

A  finite  number  of  the  diffracted  waves  represent  propagating  modes  and 
these  are  the  important  ones  far  from  the  boundary.  The  remaining  modes  are 
evanescent  and  though  these  do  not  serve  to  carry  energy  away  from  the  surface, 
they  do  play  a  vital  role  in  affecting  the  amplitudes  of  the  propagating  modes. 

The  number  of  modes  that  propagate  can  be  determined  from  the  expressions 
i  + 

for  X  given  in  Fq.  (2.30):  if  X-  is  real,  the  corresponding  mode  propsgates 
m  m 

without  attenuation,  whereas  if  X *  is  pure  imaginary,  the  mode  is  evanescent. 

in 

To  find  the  (complex)  amplitudes  of  the  diffracted  waves  we  proceed  as  follows. 

la*t  y  fCx)  bo  the  profile  of  the  surface  and,  for  convenience,  assume  the 
incident  plane  wave  to  l>e  E  polarized.  We  then  have 
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(3.5) 


(3.6) 


Having  determined  the  current  distributions  K  (x’)  ant  Kfo'),  it  is  therefore 

z  t 

a  trivial  matter  to  compute  the  amplitudes  of  ail  the  diffracted  modes,  both 
propagating  and  evanescent,  and  the  amplitudes  of  the  propagating  waves  for 
several  different  values  of  d/X,  a/X  and  6  are  given  in  Tables  IIT-1  through 
HI-4. 

A  widely  used  technique  for  estimating  the  field  which  is  scattered  by  an 
object  is  the  physical  optics  method  in  which  the  surface  field  is  approximated  by 
its  geometric?!  optics  value.  Although  the  estimates  are  deficient  as  regards 
the  polarization  characteristics  of  the  scattering,  the  method  often  provides  an 
adequate  approximation  of  the  scattering  if  all  dimensions  (including  radii  of 
curvature)  of  the  scattering  surface  are  large  co  npared  with  the  wavelength, 
and  because  of  its  great  convenience,  the  metuod  is  commonly  employed  in 
rough  surface  scattering  analyses  oven  when  all  radii  are  not  large.  A  matter 
that  is  then  ri  some  debate  is  whether  shadowing  should  be  taken  into  account, 
and,  if  s.,,  how.  According  to  the  physical  optics  method  the  current  induced 
in  the  surface  is 

K=2nAH*  (3.7) 

in  the  illuminated  region,  and 

K-0  (3.8) 

in  the  geometrical  shadow,  where  n  is  the  outward  normal  to  the  surface. 
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la  many  Instances,  however,  some  of  the  failures  in  the  resulting  estimates  of 
the  scattering  can  be  traced  to  the  discontinuity  in  the  surface  field  created  at 
the  shadow  boundaries,  and  improved  estimates  can  be  had  by  using  (3.7)  over 
the  entire  surface,  shadowed  as  well  as  lit. 

Jr  "he  present  case  of  a  periodic  surface,  the  amplitudes  of  the  diffracted 
modes  predicted  by  either  version  of  the  method  can  easily  be  obtained  by  in¬ 
serting  the  postulated  expressions  for  the  current  into  Eqs.  (3. 4)  and  (3. 6).  A 
computer  program  has  been  written  tG  carry  out  the  calculation,  and  the  ap¬ 
propriate  results  are  included  in  Tables  n-1  through  in— 4.  It  will  be  noted 
that  the  results  produced  by  either  version  of  the  method  show  very  little 
agreement  with  the  exact  data,  hi  general,  the  physical  optics  values,  with 
or  without  shadowing  included,  tend  to  be  too  small,  and  we  further  note  that 
if  shadowing  is  excluded  (or  is  not  present),  the  physical  opL.cs  esimates  are 
the  same  for  both  polarizations  (see  Appendix  A). 

It  must  be  admitted,  however ,  that  for  all  of  the  surfaces  considered, 
the  minimum  radii  of  curvature  are  too  small  to  provide  reasonable  hope  for  the 
physical  optics  method  to  be  accurate,  With  a  sinusoidal  surface,  the  minimum 
radius  of  curvature ,  pm£n  ,  occurs  at  the  peaks  and  troughs,  and  the  radius  in¬ 
creases  to  become  infinite  where  the  surface  crosses  its  mean.  For  the  four 

different  surfaces  embraced  by  Tables  UI-1  through  OT-4,  the  values  of  kp  , 

min 


are  as  follows: 

d  =  0. 2X, 

a  =  0.  IX, 

kp  .  =  0.064  , 

min 

=  1.  9A, 

=  0.25X  . 

=  2.30, 

=  0.4  X, 

=  U.2A  , 

=  0. 127  , 

=  0.2  A, 

=  0.03  X  , 

= 0.212 

To  judge  from  this  listing,  one  might  expect  the  physical  optics  method  to  be 
most  successful  for  the  second  surface  (Table  EE- 2),  but  in  practice  it  works 
best  for  the  fourth  surface,  which  uas  a  much  smaller  radius,  but  a  small 
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enough  value  of  a/X  to  be  nearly  plane. 

Before  leaving  this  discussion  of  physical  optics,  we  note  that  if  shadowing 

is  ignored  or,  because  of  the  incidence  angle  (0<  tan“*~-  ),  is  not  present,  a 

~  zira 

simple  ^Aialytical  expression  for  the  amplitudes  of  the  diffracted  modes  can  be 
obtained  by  an  extension  of  the  treatment  in  Senior  (1959).  An  outline  of  the 
procedure  is  given  in  Appendix  A,  and  the  results  obtained  from  these  formulas 

are  in  agreement  with  the  computed  values  in  the  last  line  of  each  of  Tables 

\ 

m-1  through  m-4. 

3.3  The  Current  Induced  on  a  Sinusoidal  Surface 

We  here  examine  the  nature  of  the  surface  field  distribution  for  five  different 
sinusoidal  surfaces.,  considering  first  the  results  for  E  polarization  and  then  those 
for  H  polarization,  hi  each  case,  the  modulus  and  phase  (in  degrees)  of  a  nor¬ 
malized  current  are  plotted  as  functions  of  x/X  over  a  single  period  of  the  surface 
running  from  peak  to  peak.  The  horizontal  scale  is  therefore  the  horizontal  dis¬ 
tance  in  wavelengths  and  not  the  distance  along  the  actual  surface. 

The  particular  current  normalization  chosen  is  such  that  for  E  polarization 

the  quantity  plotted  is  K  ,  the  physical  optics  value  of  which  is 

iC  z 


■/l+jf'Cx)] 


1(coS0-('(x)sine)e-Jkt,tsln8-JWcc,oe>  , 


whereas  for  H  polarization  the  quantity  plotted  is  simply  K  whose  physical  optics 

V 

approximation  is 

— jk(x  sin0  -f(x)  cos  6) 

6  6  • 


In  each  of  the  following  figures,  the  exact  computed  values  are  shown  as  circled 
points,  and  are  joined  by  a  broken  line  only  to  guide  the  eye;  the  physical  optics 
approximation  is  shown  as  a  solid  line. 
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3.3,1  E  Polarization 

For  a  surface  having  d  =  0. 2X  and  a  =  0.  IX,  the  results  for  the  three 
incidence  angles  6  =  0  (normal  incidence},  30°  and  60°  are  given  in  Figs.  3-1 
through  3-3,  respectively.  We  observe  that  for  this  relatively  small  period 
moB«  of  the  current  is  concentrated  in  the  vicinity  of  the  surface  peaks,  with  the 
current  being  almost  zero  in  the  troughs.  There  is,  indeed,  almost  an  exponen¬ 
tial  decrease  in  the  current  modulus  away  from  the  peaks,  and  the  main  effect  of 
increasing  the  incidence  angle  is  to  scale  the  curves,  leaving  the  general  shape 
unchanged.  The  phase  is  somewhat  more  sensitive  to  6,  and  whereas  the  curve 
is  almost  flat  for  6=0,  the  structure  increases  noticeably  with  6  . 

Since  kp  .  =  0. 064  for  this  surface,  it  is  not  surprising  to  find  that  the 
mm 

physical  optics  approximation  bears  no  resemblance  to  the  exact  data.  This  is 
particularly  true  of  the  modulus;  the  exact  curve  is  a  great  deal  more  simple 
than  in  the  physical  optics  one.  The  physical  optics  phase  is  also  poor  for  6  =  0, 
but  agrees  better  for  the  larger  6,  at  least  in  an  average  sense. 

The  effect  of  decreasing  the  amplitude  of  the  surface  while  keeping  the 
period  constant  is  illustrated  in  Figs.  3-4  and  3-5  for  the  case  of  normal  in¬ 
cidence.  It  is  observed  that  as  the  height  decreases,  so  does  the  current  con¬ 
centration  near  the  peaks,  but  even  for  a  as  small  as  0.01X  (Fig.  3-5)  there 
is  still  amost  a  2  :  1  variation  between  the  peak  and  trough  values.  The  phase, 
on  the  other  hand,  is  much  more  nearly  constant,  and  is  more  akin  to  that  for  a 
flat  surface  than  a  sinusoidal  one . 

Results  for  d  =  0.4X  and  a  =  0. 2X  with  6=0  and  60°  are  shown  in  Figs.  3-6 
and  3-7.  By  comparing  Fig.  3-6  with  3-1  it  is  seen  that  doubling  both  d  and 
a  has  little  effect  on  the  modulus  of  the  current,  but  has  a  marked  effect  on  the 
phase.  This  is  true  also  of  the  curves  for  oblique  incidence  (cf.  Figs.  3-7  and 
3-3). 
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surface  field  for  E  polarization  with  d  =  0.2X,  a=  0.  IX  and  t*  =  30' 
Lysical  optics. 
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.3-6:  Modified  surface  field  for  E  polarization 
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All  of  the  surfaces  so  far  considered  have  had  rather  small  values  of 
d/X  and  a/X,  and  have  not  been  such  that  one  could  expect  physical  optics  to 
provide  a  reasonable  approximation.  For  a  surface  with  a  much  larger  period, 
d  -  1. 9  X  with  a  =  0. 25  X,  at  normal  incidence,  the  appropriate  curves  are  shown 
in  Fig.  3-8.  Physical  optics  now  constitutes  a  r  .ore  accurate  approximation, 
and  though  there  are  still  noticeable  discrepancies,  qualitatively  as  well  as 
quantitatively  ,  the  physical  optics  curves  do  approximate  the  exact  data  in  a 
mean  sense.  The  exact  current  modulus  still  has  a  maximum  at  the  surface  peaks, 
but  the  minimum  now  occurs  where  the  surface  slope  changes  sign,  i.e.  where 
the  surface  crosses  the  x  axis.  The  phase  remains  fairlv  constant  over  the  lower 
part  of  the  concave  portion  of  the  surface,  but  changes  rapidly  from  positive  to 
negative  and  vice  versa  where  modulus  has  its  minimum. 

Since  the  particular  case  considered  in  Fig.  3-8  has  also  been  treated  by 
Neureuther  and  Zaki  (1969)  using  their  numerical  program,  it  is  appropriate  to 
compare  the  data.  This  is  done  in  Fig.  3-9.  The  agreement  is  very  good  in 
spite  of  the  rather  different  numerical  procedures  involved  in  the  two  methods. 

3.3.2  H  Polarization 

For  each  of  the  situations  considered  above,  we  have  also  computed  the 
surface  fields  when  the  incident  wave  is  H  polarized,  and  we  now  present  these 
data  in  the  same  sequence  as  we  did  for  E  polarization. 

The  results  for  incidence  at  9  :  0,  30°  and  60°  on  a  surface  having  d  0,2  X 
and  a  1  0.  1  X  are  shown  in  Figs.  3-10  through  3- 12.  Eight  sampling  points 
were  used  in  the  computations.  We  observe  that  the  curves  for  the  modulus 
are  roughly  sinusoidal  in  form  and  are  relatively  insensitive  to  changes  in  0, 
with  at  most  a  slight  reduction  in  amplitude  as  0  increases.  In  contrast  to  the 
behavior  for  E  polarization,  however,  the  maximum  now  occurs  in  the  surface 
trough  ,  with  a  minimum  at  the  peak.  There  is  again  no  agreement  between  the 
eJtaet  data  and  physical  optics. 

39 


with  d  -  0. 2X,  a  -  0.  IX  and  6  =  0° 
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The  effect  of  changing  the  surface  amplitude  while  keeping  the  period 
constant  is  illustrated  in  Figs.  3-13  and  3-14  for  the  case  of  normal  incidence. 
As  the  surface  height  decreases,  the  numerical  results  approach  the  physical 
optics  values,  a  fact  which  is  not  surprising  since  physical  optics  is  most  ef¬ 
fective  for  a  planar  surface.  As  in  the  case  of  E  polarization,  the  phase  curve 
is  relatively  flat  even  for  a  =  0 . 03X,  and  does  not  show  the  variation  displayed 
by  the  physical  optics  plane. 

Figures  3-15  and  3-16  show  the  results  for  d  =  0.4A  and  a  0.2X  with 
9  =  0  and  60°,  respectively.  By  comparing  Fig.  3-15  with  3-10  it  is  seen  that 
doubling  both  d  and  a  has  a  somewhat  larger  effect  on  the  current  modulus 
curve  than  was  the  case  for  E  polarization.  In  particular,  there  is  a  notable 
change  of  shape,  with  the  minimum  no  longer  occurring  at  the  surface  peak. 
Similar  changes  occur  for  6  =  60°  (cf.  Figs.  3-16  and  3-10).  Once  again,  there 
is  no  .agreement  between  the  exact  data  and  physical  optics,  either  in  modulus 
or  phase. 

The  current  distribution  for  a  surface  with  larger  period  (  d  =  1.  9  X  and 
a  =  0.25X)  at  normal  incidence  is  presented  in  Fig.  3-17.  The  modulus  has  a 
marked  oscillation,  quite  distinct  from  that  found  for  E  polarization,  with  three 
maxima  per  period.  One  of  these  is  at  the  surface  trough,  and  the  current  is 
again  a  minimum  at  the  surface  peak.  The  agreement  with  physical  optics  is 
poorer  than  for  E  polarization,  with  neither  modulus  nor  phase  being  approxi¬ 
mated  to  any  real  extent.  The  exact  data  are,  however,  in  agreement  with 
those  obtained  by  Neureuther  and  Zaki  (1960)  for  this  case,  as  evident  from 
Fig.  3-18. 

3.4  Energy  Conservation 

In  solving  problems  of  electromagnetic  scattering  from  perfectly  con¬ 
ducting  surfaces,  a  common  procedure  for  checking  the  accuracy  of  the  solution 
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FIG.  3-13:  Modified  surface  field  for  H  polarization  with  d  =  0. 2X,  a  =  0. 03Xand  6  =  0°:  — o— ,  exact; 
- ,  physical  optics. 


oduiuA 


1363-6-T 


o 


(soauiiap)  aesqa 


r)l 


FIG.  3- 18:  Modified  surface  field  for  H  polarization  wi 
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is  to  determine  the  degree  to  which  the  law  of  conservation  of  energy  is  fulfilled. 
This  is  essentially  the  check  employed  by  Petit  and  Cadilhac  (1884),  Neureuther 
and  Zaki  (1969)  and  Green  (1970),  as  well  as  by  many  others.  However,  as  pointed 
out  by  Amitay  and  Galindo  (1969),  energy  conservation  does  not  provide  a  mea¬ 
sure  of  accuracy  of  a  solution  found  by  the  Ritz  or  other  related  methods;  and 
hence,  in  order  to  use  energy  conservation  to  check  other  than  computational 
round-off  errors,  care  must  be  exercised  in  choosing  the  method  of  solving  the 
integral  equation  numerically. 

The  method  that  has  been  adopted  in  this  Report  is  a  point  matching  and 
sinusoidal  interpolation  scheme  ,  and  since  the  solution  does  not  automatically 
satisfy  energy  conservation,  we  are  able  to  use  this  as  a  check.  The  relation 
that  must  be  fulfilled  is 

N  , 

£|aJ  Re(Xm)=kcos0  (3.9) 

m=o 

where  N  is  the  number  of  propagating  modes .  The  extent  to  which  this  is 
satisfied  by  the  numerical  results  that  we  have  obtained  is  shown  in  Table  HI-5. 
For  comparison  purposes,  the  percentage  errors  implied  by  the  physical  optics 
method,  with  and  without  shadowing,  are  also  included.  Whereas  all  of  the 
exact  data  satisfies  the  energy  conservation  law  to  within  (about)  one  percent  or 
better,  the  only  case  where  the  physical  optics  solutions  does  so  is  when  d  =  0.2X 
and  a  =  0.Q1A.  Such  a  surface  is  very  close  to  being  planar. 

3. 5  Convergence  Check 

In  solving  an  integral  equation  numerically  it  is  always  desirable  to  carry 
out  a  convergence  test, 

a)  to  determine  the  number  of  sampling  points  necessary  to  achieve  the 
required  accuracy,  and 
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b)  to  test  whether  or  not  the  numerical  solution  approaches  a  stable 
value  as  the  number  of  sampling  points  is  increased. 

Although  it  has  not  been  proved  mathematically  that  the  accuracy  can  be  improved 
y  increasing  the  number  of  samples,  it  would  appear  reasonable  to  believe  so 
on  a  physical  basis. 

Figures  3-19  and  3-20  show  the  results  of  a  convergence  test  applied  to 
the  case  d=0.2X,  a=0.  IX  and  6=0  for  E  and  H  polarizations  respectively.  It 
can  be  seen  that  the  solution  does  appear  stable,  and  that  surprisingly  good 
results  are  obtained  even  with  as  few  as  4  sampling  points. 

iwo  other  aspects  of  the  computer  program  that  should  be  mentioned  are 
the  computational  time  consumed  and  the  number  of  terms  included  in  the  sum¬ 
mation  of  the  infinite  series  for  the  modified  Green's  functions.  For  E  polarization, 
a  typical  figure  for  the  computational  time  on  an  IBM  360  machine  using  12 
sampling  points  is  about  60  seconds.  This  time  includes  the  computation  of  the 
amplitudes  of  the  diffracted  modes,  the  energy  conservation  check  and  the  cal¬ 
culation  of  the  physical  optics  results.  A  typical  figure  for  H  polarization  using 
8  sampling  points  is  28  seconds .  The  time  increases  roughly  as  the  square  of 
the  number  of  points,  and  most  of  it  is  eaten  up  in  the  computation  of  the  matrix 
elements. 

Except  for  the  large  period  surface  (d=l.  9X,  a=0.25X),  where  the  number 
of  terms  retained  in  the  Green's  function  series  was  6,  only  the  first  three  terms 
were  required.  Tests  that  were  performed  indicated  that  three  was  sufficient. 
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FIG.  3-19:  Modified  surfaoe  field  for  E  polarization  with  d  =  0.2X,  a  =  0.  IX 
9  =  0°  computed  using  4,  10  and  12  sampling  points. 
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IV 

CONCLUSIONS 

We  have  here  considered  the  problem  of  a  plane  electromagnetic  wave 
incident  on  a  perfectly  conducting,  two  dimensional  periodic  surface,  and  have 
developed  numerical  procedures  for  the  direct  digital  solution  of  the  integral 
equations  for  the  surface  fields.  By  using  special  summation  techniques  fol¬ 
lowed  by  the  subtraction  of  the  dc  term  to  improve  the  convergence  of  the 
series  for  the  modified  Green's  function,  a  relatively  efficient  procedure  has 
been  arrived  at,  and  this  has  been  programmed  for  a  computer. 

Data  have  been  obtained  for  the  surface  fields  on  several  different  sinusoidal 
surfaces  when  illuminated  by  E  and  H  polarized  plane  waves  at  oblique  as  well  as 
normal  incidence.  It  is  found  that  the  polarization  has  a  marked  effect  on  the 
field.  Since  most  of  the  surfaces  considered  had  relatively  small  values  of 
d/X  and  a/X,  there  is  little  agreement  with  the  physical  optics  estimate  of  the 
surface  field.  However,  this  is  also  true  for  the  one  surface  of  much  larger 
period  that  was  examined. 

Knowing  the  surface  field,  the  amplitudes  of  the  diffracted  waves  in  the 
discrete  angular  spectrum  representation  of  the  scattered  field  can  be  computed, 
and  this  has  been  done  using  the  exact  surface  fields  as  well  as  the  physical  optics 
estimates  with  shadowing  either  included  or  Ignored.  Here  again  the  physical 
optics  predictions  are  deficient,  and  whereas  the  results  derived  from  the 
numerical  program  satisfy  the  conservation  of  energy  law,  the  physical  optics 
values  do  not. 

In  the  continuation  of  this  study,  the  numerical  program  will  be  applied  to 
a  wider  variety  of  periodic  surfaces  in  order  to  build  up  an  understanding  of  the 
surface  field  behavior,  including  its  dependence  on  profile  size  and  shape,  and 
the  direction  and  polarization  of  the  incident  field.  It  is  our  hope  that  this  will 
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enable  us  to  develop  synthesis  procedures  for  the  efficient  computation  of  the  fields 
scattered  by  periodic  surfaces  of  large  amplitude  and  period,  leading  ultimately 
to  better  prediction  techniques  for  rough  surfaces  in  general. 
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APPENDIX  A:  PHYSICAL  OPTICS  APPROXIMATION 

Although  the  main  purpose  ot  this  report  is  the  development  of  a  numerical 
teclmique  for  the  determination  of  exact  values  for  the  surface  and,  hence, 
scattered  fields,  we  have  found  it  desirable  to  compare  the  results  obtained 
with  those  provided  by  the  physical  optics  approximation.  This  approximation  is 
a  physically-based  one  which  postulates  an  explicit  form  for  the  surface  field 
arrived  at  by  assuming  each  element  of  the  surface  to  bear  that  current  which  it 
would  were  it  part  of  the  local  tangent  plane.  The  calculation  of  the  scattered 
field  is  then  reduced  to  quadratures. 

In  many  instances,  however,  and  a  periodic  sheet  is  one,  an  analytical 
evaluation  of  the  physical  optics  integral  is  a  difficult  procedure,  particularly 
in  such  cases  where  part  of  the  surface  is  shadowed,  and  some  of  the  shortcomings 
of  physical  optics  estimates  in  general  certainly  arise  from  sloppiness  in  the 
evaluation  of  the  integral.  In  still  other  cases,  the  physical  optics  result  proves 
to  be  more  accurate  if  shadowing  is  ignored*  (see,  for  example,  Adachi,  1965) 
and  it  is  of  interest  to  observe  that  for  a  sinusoidal  sheet  it  is  then  possible  to 
produce  an  exact  evaluation  of  the  integral.  Needless  to  say,  however,  the 
resulting  scattered  field  is  still  subject  to  the  unknown  errors  inherent  in  the 
use  of  the  physical  optics  approximation,  and  to  the  neglect  (if  present)  of  all 
shadowing  effects . 

The  procedure  is  directly  analogous  to  that  given  by  Senior  (1959)  for  the 
particular  case  of  an  H  polarized  plane  wave  at  normal  incidence  on  a  sinusoidal 
sheet  and  consists  of  three  steps: 


* 

This  modified  method  is  acmetimes  called  'extended  physical  optics'  . 
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<i)  writing  down  the  physical  optics  integral  for  the  scattered  field  without 
any  shadowing  included; 

(ii)  asymptotically  evaluating  this  expression  appropriate  to  an  observation 
point  at  large  distances  from  the  sheet;  and 

(iii)  matching  this  expression  to  a  discrete  angular  spectrum  of  waves  to 
obtain  their  amplitudes.  Since  these  amplitudes  are  independent  of  the  field 
point,  we  have,  in  effect,  produced  an  exact  evaluation  of  the  integral  valid 
certainly  in  the  half  space  above  the  sheet. 

Let  us  take  the  equation  of  the  perfectly  conducting  sheet  to  be 


y  =  a  cos  Kz  (A.  1) 

where  a  is  the  amplitude  of  the  corrugations  and  2jr/K  =  d  is  the  period.  If 
the  incident  field  is  E-polarized,  we  write 


i  a  -jk  (x  sin  6  -  y  cos  0  ) 
E  =  z  e 

(cf.  Eq.  2.20),  implying 


(A.  2) 


„  A  A 

-  -  Y  (cos  8x+  sin  8  y 


-jk  (x  sin0  -y  cos  6) 
e 


(A.  3) 


where  Y  is  the  intrinsic  admittance  of  free  space  . 

By  virtue  of  the  periodicity  of  the  sheet  and,  hence,  of  the  problem  as  a 
whole,  the  scattered  field  can  be  expressed  as  a  discrete  spectrum  of  waves 
which  waves  are  certainly  outgoing  as  regards  y  >  a.  Thus,  we  have 

oo  -jk(x  sinfl  +  y  cos  0  ) 

ES=£  £  Ame  m  (A.  4) 

m=-oo 

where 
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k  sin  0  =  m  K  +  k  sin  0  , 

m  _  1 

k  cos  0  ={k^  -(mK.  a0)^  , 

m 


(A.  5) 


and  that  branch  of  the  square  root  is  chosen  having  imaginary  part  non-positive. 
Application  of  the  physical  optics  approximation,  followed  by  steps  (!)  through 
(ill)  above  now  gives 


,  .m  k  la  mKa  .  .  ,  .  . 

A  =  -j  _ eosfl-" —  sin  0  J  (u  ) 

m  _  ^  m  m 

Vk  -(mK-Hcsin0)  ' 


where 


/  f~2  2*\ 

H  -■  a  (jc  cos  0  +V  k  -(mK-tfcsin  0)  ) 


In  particular. 


A  =  -  J  (2a  k  cos  0  )  , 
o  o 

and  as  a  0  , 

A  1,  A  — ►-  0.  m  =£0  , 
o  m 

in  agreement  with  the  known  solution  for  a  flat  sheet. 
If  the  incident  field  is  H-polarized,  we  take 


(A.  6) 


(A.  7) 


i  *  -jk(x  si»0  -v  cos  0> 

H  =  z  e 

implying 

i  1  .  a  .  .  _  a.  -jk(x  sin  0  -y  cos  0) 

13  =  —  (cos  0  x  +  sin  0  y)  e 

and  expand  the  scattered  magnetic  field  as 


(A.  8) 


(A.  9) 


s  s 

H  =  z 


UJ 

I 


-jk(x  sin  0  +  y  cos  0  ) 

m  m 


m=-eo 


(A.  10) 
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whore  9  is  as  before.  On  evaluation  of  the  physical  optics  integral,  we 
m 

now  obtain 


A'  =A 


(A.  11) 


m  m 

The  equivalence  of  the  results  for  the  two  polarizations  is  consistent  with  the 
known  fact  that  the  physical  optics  approximation  for  a  perfect  conductor  is 
inherently  polarization  insensitive. 

Results  computed  from  Eqs.  (A. 6)  and  (A.  11)  agree  with  those  given 
in  the  appropriate  columns  of  Tables  m-1  through  EI-4  and  obtained  from  a 
numerical  evaluation  of  the  integrals  in  Eqs.  (3.4)  and  (3.6)  . 
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APPENDIX  B:  COMPUTER  PROGRAMS 

Two  separate  computer  programs  have  been  written,  one  for  E  polarization 
and  the  other  for  H  polarization.  Each  program  consists  of  the  following  parts. 

1.  Main  Program 

The  functions  which  this  performs  are 

a)  reading  in  the  input  data  consisting  of  the  amplitude  and  period  of  the  surface 
(measured  in  wavelengths),  the  angle  of  incidence,  the  number  of  sampling 
points  and  the  number  of  terms  employed  in  the  series  summation; 

b)  calling  the  subroutine  which  computes  the  matrix  coefficients,  including 
the  analytical  treatment  of  the  singular  term; 

c)  computing  the  physical  optics  current; 

d)  calling  the  subroutine  which  inverts  the  matrix  and  hence  obtains  the  induced 
surface  field;  and 

e)  printing  out  the  appropriate  numerical  results. 

2.  Subroutines 

a)  Subroutine  CAVE 

This  determines  the  appropriate  constants  in  the  sinusoidal  interpolation 
formula  and  calls  the  integration  subroutine  CSIMEQ  to  integrate  the  kernel 
function  over  the  sampling  cell. 
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b)  Subroutine  FUN 

The  constants  in  the  sinusoidal  interpolation  formula  are  here  computed. 

c)  Subroutine  SUBC 

This  computes  the  kernel  function  of  the  integral  equation. 

d)  Subroutine  CIMPS 

The  numerical  integration  is  performed  using  the  Newton-Cotes  method. 

e)  Subroutine  CSIMEQ 

This  performs  the  matrix  inversion  using  a  Gauss  reduction  method. 

f)  Subroutine  ENERGY 

The  amplitudes  and  phases  of  the  diffracted  waves  are  computed  and  the 
scattered  energy  found. 


67 


I 


I 

r 


(2 

!o 

N 

II 

pc 

c 


ll 

£ 

!± 

1C 


p 

I 

s 

II 

ft 

c 


1 

>  1 

ml  1 

k-  ' 

to 

O  I 

o  « 

it  ! 

O  *s. 

it  — 

»  4- 

oe  o 

2  < 

o 

i 

It  CM 

•  /'.  k- 

“J 

a  — 

o  <t 

2 

2 

-  o 

O 

<  M 

to  >. 

u 

'  tm 

40  X 

► 

X  1  «> 

►  z 

oe 

o 

to  3 

O 

a  O 

(J  It 

k- 

a.  cm 

•  » 

to 

.  w 

•M  > 

•> 

tr  lu 

-*  O 

CD 

2  jlO 

a  o 

< 

t-  < 

•  -v 

o£  ;X 

—  m 

uj  !a 

a  < 

o 

k-  * 

•  K 

a 

I- 

X  < 

► 

<  C 

X  o 

H 

O  CM 

*  >» 

* 

z 

to  —  — 

•“  Ui 

*-t 

o 

c  o 

X  X 

X 

*-« 

X  UJ  <  CM 

< 

r 

H 

it  >  X  w 

• 

to 

< 

<Nk 

CM  X 

» 

K 

u.  2 '  •  40 

<  • 

Cl 

r 

C  -  C 

k-  *-> 

H 

a 

UL'OC 

<  » 

X 

>  fM  • 

a  < 

X 

o 

C  111  - - 

*M  U. 

2 

o 

C  I  M  O 

z 

3 

k  (-  *  CM 

<  * 

a 

oe 

CO  •—  — 

k-  C 

• 

c 

at  c  <i 

a  to 

e; 

u. 

a  c  <v  X 

X  a 

< 

I  U.  -  N 

k-  ► 

Ci 

IO 

k*  2  • 

♦  c  •• 

z 

k- 

U.  ►  — 

X  k- 

> 

z 

X'  z 

*  o  z  c 

x  a 

C!  X 

» 

< 

C  <  C  (M 

•  » 

u 

* 

»  z 

k- 

u.  a  i* . 

a  u' 

• 

z 

» 

IO 

t*  *■'  to 

X  2 

a 

► 

c  < 

z 

x  =  a  <  <  » 

X 

< 

k-  k- 

o 

to  •  x 

►  i/: 

it 

H 

U. 

u 

oe  *»  a.  z  * 

u 

u.  C  X 

CO 

c 
<c , 


t-  I 

I 

*  be 

*  C 
-»  it 
io  t- 

< 

S; 

*  5 


C  O  c  •  • 

C  *—  IV  c 

Of  c  •  C  x 

a  c  c  ev  v 

■»  <\  —  ♦*■ 

z  a  -  c  < 

•-  UJ  «  a.  I- 

<  a  *,« 

S  »2  -  c 

C  C  V 

U)  c  *" 

XZie-'Z 
I-  ~  2  I-  C 

k  UU  «/■  X 

f  c  *  *-  5 

U.  C  *-  c  c 

_■  c  c  •  c 

*-  «/■  2  — 

a  *-  c 

22C 

e  c 

*  o 


w 

to 


I 

k- 

io 


I-  IT 

•  k» 

_  a  (v 

a  -c  x  • 
z  c  <  ir 
•  z  <  *■»  *  • 
a  p  •'  oc 


u. 

>- 

» 

2 


zlc 
k-  I 


u 

IT1 

to 

c 

•  2 

» 

r-  (V 


♦  ecu 

<C  X  ]K 

->  e 

x  «•  *  • 

UjZr. 

•j  of  !z 
a  aio  » 

*  »-  «| 

C  x  a  it 

o  it  a  !*- 


x  c. 


it 
■  fO 


Z  z 

«-  oe 
a  c 
a.  u. 


UJiO 
k"  1  • 

«'c 

Jl  It 

r> 

or 

j,o 

<  « 

o 


iB6348-T 


in 

*e  i 
m  ! 
<r  ] 

IT  k-j 
«-  a 

*•*  •: 
•  (M 

rr  M 

H 

a  a 


it  it 

c  c 
,#  « 
!*»•  ** 
2  »  •*» 
xj&  a 
v  —  — 
C  V‘  2 
X  C  k- 
Z  n  O  to 
b  _j  p  w 
Z  a  ^  / 
x  ct-io 


a  a 

«  * 

n  *  to 

*-  *■»  to 

U)UI  « 

I  X  z 
•  i-  i-  to 

I/!  ♦“  C  • 
O  to  O  rv 

«  n  Jl  « 

|e  to  c  o 


68 


IM  V 

f  -* 

•  It 

*-  e 

*-  M 

a  < 

H  H 
O  ml 
to  a 

a  c 


i 


136^-6-t 


* 

M 

ft 


I/) 

+ 


rv 


X 

1 

« 

< 

a r 

1 

5 

i- 

UJ 

-j 

X* 

t 

1 

UL 

I 

c 

> 

I 

— » 

1 

k- 

! 

Iftk 

• 

ft 

•“« 

•  < 

^4 

< 

or 

♦ 

• 

» 

k- 

< 

1 

» 

• 

ft 

-J 

-i 

• 

o 

ft 

O 

U. 

r 

o 

ft 

<— 

c 

tf 

i 

*»» , 

# 

o 

* 

— ’ 

z 

# 

w 

*— 

c 

IH 

ft— • 

ft— 

ft 

c 

1/) 

| 

ft— 

*— 

ft— 

•J 

a 

a 

►— 

a 

< 

UJ 

v. 

V 

a 

V 

ft 

X 

►— 

»— 

'v 

>— 

<1 

k- 

X 

X 

X 

1- 

V 

Nv 

X 

-1 

z 

e 

c 

*v 

C 

U 

*-* 

i 

X 

X 

IV 

c 

X 

c 

1 

+ 

1 

C 

X 

1 

ft 

z 

r 

o 

rv 

4 

• 

C  I 

cv 

<v 

— 

ft- 

4 

IN! 

<  IV  C 

4 

< 

K  l 

ft 

c 

1 

1  ft 

K 

< 

1 

ft 

« 

ft 

k 

ft-* 

IV.  ft 

1 

00 

X 

<  • 

< 

c 

0- 

< 

9 

X 

C 

%*» 

w  #■» 

c 

< 

w- 

0 

*■* 

X  C 

1- 

1-  X 

w 

K 

• 

or 

k 

a 

or  ft- 

A 

1- 

or 

• 

IT 

Q 

C  vr 

or 

o 

c 

<s. 

4 

ir 

or  4 

<4 

o 

«/. 

w 

UJ 

c 

«V 

-v 

v.  C 

< 

IS 

v. 

ft 

*■* 

X 

< 

• 

•  X 

V. 

• 

0- 

or  •* 

V 

«-* 

0,-  V 

« 

• 

9mm 

UJ  • 

Z 

►— 

• 

— *  ft— 

a1 

— 

a 

00  c 

X 

a.> 

4 

1  X 

-i 

4 

1 

V 

0 

** 

ft 

i- in 

—  ft 

• 

— 

ft- 

C 

Uu  •  X  II 
i  e  X  *-  . 
*■*  —  •» 
IlC*-4 
H  ^  |« 

<  r  u  •» 
u  r  k  c  >• 
or  </  x  C  x 


C 

k- 


O  a 
ts  I 
a  w 

II  H 

rv 

<  < 


•  X  C  X  *- 
25  IV  —  *-* 
*“  1/  1/  IX 

<  «  c  II 


x 

u-rc 

*-  or  e 


X  II 

r  * 

or  < 

0-  <*■ 

c  c 

i\  rv 


«  r-  «.  iu  ^ 

—  rv.  “  r  ♦ 

«  v  v  x  > 

II  c  II  —  k- 

II 


W 

U 

c 


X 

UJ 

X 


I  i 


X 

UJ 

c 

z 


a 


X  X  Z  C  «J 

:c:c<u 
v  e  «/■  c  c  e 

IV 

c  rv 
rv 


» 

C 

z  z  • 
•  •  c 


it  ii  ii  r  h 

—  X  *—  Z  OL  >— 

x  r 

^  IV  •  H  |T 

_  2  -J  It 

c  c  -  c  _j  c  ►. 

C  C  c  o.  C  X 

u. 

rv 

c 


c 

j*. 

b~ 

la 

\n 

X 

0 

,< 

z 

• 

ft 

IV 

,X 

c 

V 

1- 

0 

o 

1- 

a  or 

« 

z 

C  C 

o 

0m 

k-  k* 

z 

0 

1/  00 

tM 

>  it 

or 

-  C 

U* 

ft 

e  iv 

ff  o: 

< 

• 

i r. 

•c 

a  C 

r- 

a' 

—  k- 

4  4 

— » 

•  | 

c 

<— ' 

or 

IV.  ■  » 

*0 

x  c 

c 

\  »■  • 

or 

•J 

*  c 

i  1 

K 

—  —  c 

K 

u. 

X 

<—  k- 

tr 

►*  X  — 

< 

c 

*—  00 

■—  w 

—  ft  — 

X 

ft 

0-  X  0 

k-  H 

IT 

k-  c  a 

— 

X  mJ  0- 

U-  V 

er 

or  —  x 

UJ  Z 

•J  IT 

ft  uj  or  • 

c  c 

c  a  uj  z 

X  • 

a  • 

•  c  a  c 

c  c 

<  c 

c  —  c  ► 

♦  I 


c  x 
«c  |  •-* 
—  m  ^ 

X  X  II  x 

II  -  -  I,  _• 

•“  H  —  X  _ 

x*  4] 

x  X  K  X  u 


Xu.ll 


•w  k- 
it 
a  *- 
-  C 


n  C 

—  IV 

-  C 
k-  ►“ 
it 

C  C 
c  c 


IT 

c 

c 

IT 

c 

IV 


IL 

c 

ft 

m* 

IT 

H  V.  ft  »»■  C 

““  C  IM  II  I 

C  X  X  X 

—  H  II  X 

*"  0,‘  —  1*1  X  — 

^  Z  *“  H  H 

~  1 - C  x  * 

Ctricxx 

< 
c 
rv 


c 


60 


I 


I* 

:x 

* 

o 


1363-0-T 


i  • 


c 

< 


z 

| 

;  : 

o 

m* 

m* 

*— l 

1 

| 

i« 

■  » 

. 

<✓) 

! 

! 

o 

|X 

# 

| 

* 

1 

1 

UJ 

1 

; 

'ft. 

X 

1 

UJ 

1  ► 

,4 

1 

H 

‘ml 

| 

w 

; 

;4 

« 

i 

| 

a 

UJ 

ft 

X 

,« 

ft  ft 

(Si 

j 

-  UJ 

ft  ft 

# 

t 

1 

o 

w 

ft  ft* 

« 

1 

* 

; 

, 

,x 

r~  «" 

ft 

J 

— 

!«C 

o 

mm 

l 

ftl 

:  • 

mm 

3  a 

mm  1C 

1 

« 

X 

X  X 

« 

a 

. 

ml  J 

i-  * 

! 

• 

N 

4  « 

u  a 

! 

•» 

UJ 

ft. 

a  a 

*  — 

a— « 

UJ 

It 

a  af 

—  a 

j 

X 

Of 

» 

N  ^ 

x  — 

« 

c 

X 

ft  ft 

x  z 

! 

c 

UJ 

N 

ft  ft 

#  — 

j 

o 

» 

ft  ft 

c  un 

I 

& 

2 

IT 

ft  ft 

i-  * 

«■»  1 

— 

C 

i  . 

. 

3  C 

a  u 

m\ 

z 

i/j 

c 

ft  a 

-  X 

♦  • 

• 

• 

►a 

ft- 

r- 

o 

ftj 

e  - 

It  4 

M 

<v 

«r. 

z 

J 

a 

—  a 

<  c 

C  - 

•  ’ 

» 

• 

ft 

ml 

“*■ 

_j 

2  • 

x  < 

c  + 

px ' 

ft 

a 

c  - 

vn 

a 

*  X 

ft  X 

*  . 

» 

4ft 

• 

X 

^  3 

a 

ft*  fft 

"1  - 

a  r- 

H  ‘ 

H 

>• 

4 

ftj 

a 

xj 

M  . 

ft  4 

x  ~ 

2 

2 

z 

mm 

t-  M 

'X 

ft.  (Si 

*  r 

<  — 

*4 

*•* 

*— < 

♦ 

a  • 

V 

IA 

•  * 

4  Z 

14  4  4 

* ; 

» 

* 

« 

O  OC 

a 

ft. 

— 

*-  < 

X  O  4 

►  > 

> 

»■ 

m  ' 

• 

4 

1  • 

*•  a 

4  *- 

hUu'C 

“ 

C  3 

vr 

X 

w  » 

0  4 

Vt  V  W  »- 

c 

c 

►- 

►"  * 

4 

3  X 

—  ft  # 

#  —  a 

<J  . 

< 

4' 

a 

4 

u 

» 

•  C 

ft  ft  ft 

x  a  x  C 

• 

» 

» 

a 

X  • 

•* 

U 

ft* 

-*  ft- 

ft  &  ft 

x  a  C 

<\ 

<s 

ft 

*/. 

*“  2 

a 

ft-  • 

—  v  a 

—  •  <j 

a 

#r 

4 

♦ 

a  • 

2 

2 

-  r 

3  C  * 

«  <  «  — 

•  • 

• 

» 

m~ 

H*  Z 

ft 

•sj  » 

a  c  • 

w 

►  -  <  ft 
I  •  •  ft 

a.  *-  c  a 
I  C  v 

UkJ  •  ft  — 

►■.II 
4  c  *  ►“ 
—  c 

—  H  *  I 
C  4  -  I 

-  «  C  H 
<  eft 

c  ~* 


•  ml 

O  UJ 
a  c 

•  I 

l 

x  x 

—  x 
« 

U.  ft 

—  a 


CC*  4  ff 

ml  —  «■(>-»* 

a  or  Ui  — »uj 
C  >  C  UJ  > 
X  4  <«  X  C  4 
♦  I  I  <J 
X  C  X  X 

X  «J  ►*  X  X  ml 

»  ml.  *  *  ml 

^  <  C  ^  ft  <1 
ffeettfu 


JUJU.' 

Uj  c  >  ft 
C  X  *  X 

♦  ♦  c  « 
V  X  C 

X  X  ml  ►** 

•  II  ft'  0 

•“  4  a 

*  a  v  *- 


-  «  -2  3 

*  i  o  *  — 

•—  UJ  W  c 

—  U  X  »  -J 

«  UJ  I  ► 

ft  3  »-  w 

~  *  .  c  **• 

ft  ft  ►“  IT 

C  Jft 

UJ  -J  • 

>  4  C 

2  U  C 


at 
►  X 
C  K 
cr 

Z  Uj  Uj  l- 

CS2i 

C’ZujC 
n  «“ 

-  ^  -j»- 
-*2 
ft  C  4  •- 

3  e  c  v 
a 

ft 

»r  C 
<\  H 


X 

c 

ft 

ft  — 
t*' 


—  X 

-  c 
» 

ft  W 

fft 


*  *  ~ 

*-  jE  ► 
or  C  oc 
a  a  « 


m 


2  X 
•*  or 
Or  C 

a  a 

|t 

<*• 


^  w 

*  cc  v> 

*  4  * 

•*e  * 
«  M  c 

ft*  ^ 
w  ft* 

c  *** 

*  41 

C3I 
C  K  ft 


4  C 

ft  4 

ft  ft 
M  ft 

—  «  a  < 

rs!*" 

£  it  £  £ 

4  4**“ 

X  I  c  * 
a  ft  C  ft 

»r 

it 


70 


1M3-6-T 


00 

O  V 

X  4 

00  K 

o  < 

•  c 

*/>  *v 
OO  K 

*  2 
IT  ~ 

U  u. 
•  » 
—  > 

—  C 

Q.  u 

*  S 

—  m 

a  < 

*  i- 
*  4 

X  C 

*  S 


1 


4 

l«4 

^  ! 

N.  • 

K 

X 

K  . 

N  X 

ai 

Ui 

Hi  ! 

—  4  » 

O 

o 

o 

U  h*  ►- 

2  <  • 

z 

X 

X 

*  C  4 

• 

» 

• 

tr  -v  a. 

«P* 

ff* 

2  2 

4 

4 

4  , 

♦  <  » 

• 

>- 

<  ►-  C 

t r 

• 

• 

► 

2  u.ir 

< 

2 

2 

2 

*  la 

2  ¥-  . 

• 

a 

a 

o  . 

a 

*"•  •  C 

IN.' 

♦ 

• 

•  ! 

♦  a  — 

< 

*• 

I- 

v-  —  a  a 

2 

2 

2 

<_  c  •  * 

•M 

•  ^  a  v. 

**■» 

• 

• 

» 

*  —  *  2 

4| 

► 

>■  »-  ft,  fv  «rt  *r 

Ci  £  < 

>* 

—  <41  < 

* 

V.' 

V. 

C  C  C  f  X  ►» 

<  —  2  • 

X 

* 

4 

4  ♦ 

<  c  *  o’ 

2 

• 

• 

•  IV  ft  —  «v 

—  f\  C  o' 

ft 

ft 

fv  4  C  c  «-  *-  44  -  4 

u. 

a 

<t 

*  *-  C  C  ♦  1  »-  »-  »  *»■ 

>  C  **  * 

» 

• 

• 

•  4  ¥  X  44  X  4 

<  f\  J 

4* 

4* 

^  ^  —  •  ft  ^  •  4*. 

t_-  W  4  K- 

4  ►»  tr 

5? 

« 

C 

C  If;*-  2  -  -  M  -  «. 

u.  <  * 

»  u. 

tr 

00 

*f*  *»  •  •  4  4  i-  —  4  t~ 

2  2  t*  2 

«C 

a 

a 

ft  —  O  0  »  K  <o-  —  1  . 

—  c  »- 

-J 

X 

X 

1  4  Hi  lo  <•  -  44  —4 

X  4 

VJ 

*■*» 

•*> 

*“  I  *  «  fft  B  1  (f  4  ■ 

— •  \s  2  • 

u.  2 

U  2 

4  C  * 

C  C.' 

W  —  4  |  —  —  |  _ 

c  2  c  oo 

e  4  K 

II 

T»  *  X 

X 

KfUffCMCI* 

#  U.  I  o 

a  u  4  * 

> 

» 

B  W  • 

-J  » 

» #  **-w  *  *  ►-  »  • 

'lie 

I  ►“  x  ■ 

4  «»  4 

4 

•J  —  —  «—  «—  *•—  •—  •_ 

c  ♦  J 

C  x  _  ft  u  * 

U  4  u. 

4  U. 

*T  C  t-  • 

^  A 

c  «■  c  * 

4 

*— 

2  <-  2 

C  X  L  <  •»  >»  4  «  C  «  «  »J  B  * 

72 


1363-6-T 


"  1  ‘T  ■  ■'-** 


IS. 

o  4 

Q  "v 


•  4 
2  4 

IS  ft 

o  < 

►  c 

IS  ^ 
i/"  X 

•  2 
^  C 
C  a 


—  > 

—  C 

a  o 

•  \ 

—  rr 
a  < 

•  »- 
*  « 
x  C 
•  ■v 
->  u. 
X  x 
< 
V  * 
fV  x 
4  » 
ft  »- 
4  • 
C  4 
•v  ft 
2 
4  » 
*-  C 
u  >s 
I  a 


•  c 

i  — 

—  >  a 

-  c  •  * 

X  fv  ft  l_- 

-  -  *  I 

2  *  «  - 

~  »  •  i/ 

a.  2  * 
C  *  vr 
2  <\  C  «/ 
C  •  x  ft 
•»  v  • 
ft  is  —  2 
C  ~ 

2  4  ft  v' 
2  4  • 

ft  2  C  I 
c  V  K 
X  —  V-’ 

ft  ^  2  • 

_■  2  C  *r 
Q  U.  I  \S 

size 
c  ►-  c  . 

V  C  «- 


c  c  c 

ft  >-  ft  ft  X 

c  * 

c  c  c  c  c 

•  C  c  C  ’*• 

X  **■ 

z  -  -  -  r 

r  -  cv  x  2 

u.  -  X 

•  te  *  •  ♦  — 

x  x  e  c  c  £ 

•  i  u.  ft  ft  C 

«  x  •  »  .  ft' 

X  «  4  4  4  */■ 

ft  **■  u  ft  ft  ft‘ 

--2  22  - 

a  a  —  w  -  k 

X  *  ~ 

2  C  3  ft  ft  ft  2 

•  X  —  ►  »-  ft 

rv 


i/>  i/i 

*/>  «/) 

tt  * 

•  • 

<M  <VJ 

•  * 

2  2 

i/i  *r. 

c  o 


2 

x 

IS 

C 

c 

«/: 

✓ 

4 

2 

5 

ft' 


2 

x 

%S 

c 

ft' 

« 

%S 

\S 

4 

2 

IS. 

ft' 

« 


X  X 

2  2 

—  X 

>  2 

</  *»' 

-  r 

4  I 

v  u. 

V  *►  ✓ 

C  ft 

(T  f  V 

Z  »» 

I  1 

*  2  2 

C  -  C  ~  ft  s 

ft  x  ft  *  —  n  T 

»-  •*  x  —  2  x  C 

C  2  C  2  *4Zft2 

CftCft>ftft°ft 

rv  <*■  «» 


73 


1363-q-T 


c 

> 


x 

< 

X 


z 

X 


3 

C 

a 

z 

X 

II 


c 

*£  UL 

X 

C  ►- 

X 

c 
C  X 
e  o 

»  3 
C 
X  Of 

.  X 

C  X  z 
a  ii  * 

.  x  * 

ii 

_  u_  — 

-  X 

—  *-  ir 

a  X  C 

*-  c  c 


3 

If 


c 

o 


4/) 


«-  c  *- 

•  <  II  *-> 

•  II  —  — 

’  — .  •  if  II 

•  •  *-*-*- 

!  3*  •—  cc  ac 

—  II  ~ 

<  <  |_  o 


z 

X 


C  c 
X  c 

C  u. 

X  X 
X 

c 

c  u. 

c  - 

—  if 

z  - 

a  C  X 

—  Of 

•  x  —  a 

o  < 

a  u  ”5  ii 

•  C  —  - 

C  X  - 

II  “J  Of  — 

-  u.  a  w 

—  X  ”5  X 

—•  H*  II  Q( 

_  a  *-  a 

a  -  a  - - . 

a 

C  c 

r- 

\S 

c 


z 

X 

z 

X 

X 

c 

Ml 

II 

_J 

X 

X 

c 

4** 

* 

c 

U.' 

”9 

X 

X 

* 

►1 

X 

X 

X 

X 

X 

w 

1/ 

II 

z 

cc 

•m* 

*  - 

" “ 

c  -  X 

«£ 

u. 

\S 

X  X  — 

fV 

x 

c 

*  a. 

X 

* 

X 

Of 

c  ►*  # 

c 

# 

a 

e  —  — 

< 

c 

M 

«  X 

a 

• 

w 

z 

*»• 

-  i  .  - 

<4 

< 

< 

a 

z  —  • 

X 

«v 

C 

z  z  -  -  c 

ec 

»  m 

X 

• 

•  »  *  *  <  • 

^  i— 

C 

X  z 

fr— 

Z 

O  C  C  1  • 

enz 

•  Z 

a  — 

II  * 

C  X 

u.  a  a  —  —  c 

n  » 

C  1 

y  x 

II  X 


II  if  *”•  — 

-  II  if 

if  a.  i-  * 

-  X  *- 

X  I-  <  - 

a  < 

6.  x  C  II  -  —  X 

"V  c  c  c  «  <  V 


n  C 
*5  e 
»  ii 

>-  a 


III  Xr 
X  Ii  «■* 
Z  uj  *-  ♦ 

?  *r 


3 

tf 


<4  x 
a 

c 


x  x 

li  n  ii  it 

—  ”s  — ■  ec  —  t 

n  ii  a:  ►»  x  — 

X  CC  *  —  •  X  X  g-  | 

—  Z 

aCC - CZC  H 

*-cc<)cc«ic:'-cx 

< 

c  Of  x  i-  x 

x  ec 

C 


v"  Z 

I  • 

—  if 

li  +  lb  if  Uj  it  li 

«.  .C“5XXwXx  — 

♦  C  a  X  Z  <r  Z 
X  —  X  C  if  •-  li  —  a:  •-* 

||  II  w  r-  II  X  —  X  I  r- 

C  X  X  Z  X  Z  X 

aXuCrCwC  c 

yt/t-Ci/’v.'xcaC 

u. 

*  c  <  <T  C 

IV.  —  (\  a 

c 

c 


N  II  W 

-  X 

x  a 

ii  cv.  a 

-  ii 

-  c  _ 

a  C  - 


or 

••  z 

—  a 

»-  X  C 

—  u.  Z 

x  a  a 


77 


«FA|  FUNCTION  CA«S1(7) 


Z 

u. 

a 

x 

<_■ 

if 

c 


136^-6-T 


c 

< 

•— 

u 

If 

z 

> 

a 

** 

I 

X 

^  p"" 

a 

a.1 

— « 

«— « 

ii  cc 

X 

u 

<v  * 

a 

X  • 

a. 

K 

* 

*  p~ 

*  IT 

*  pp. 

•1 

**• 

r- 

* 

<v 

X 

< 

X  UL 

p- 

ft 

IL 

w  • 

CL 

X 

tt 

X  fs. 

tt 

U_  it 

»- 

> 

—  a 

If  < 

< 

It 

f\  W 

<i  t/> 

— 1 

c 

x 

a  c 

Ot  — 

X 

—  <  — 

•  u 

< 

c 


a  it 


—’  2  <\  —  — 


< 

a 

c 

C 

a 


a 

X  Z 


w 

c 

< 

-tt 

X 

IT 

a 

< 

* 

a 

• 

c 

a 

X 

a 

4 

c 

IT 

<w' 

w 

X 

♦ 

<\ 

X 

c 

< 

•»> 

« 

■tt 

2 

II 

■tt 

—• 

< 

2 

CL 

Q 

c 

« 

2 

< 

rv 

5 

i/ 

p— 

* 

cr 

• 

» 

• 

< 

tt 

a 

4 

< 

— - 

IT 

t 

t 

a 

X 

fv 

w 

X 

• 

c 

2 

2 

X 

tt 

2 

IT 

< 

*- 

tt 

< 

<\ 

r* 

sr 

ft 

a 

H 

a 

u. 

u. 

•— 

c 

tt 

♦ 

tt 

a 

X 

< 

• 

<\ 

c 

X 

X 

X 

• 

CL 

w 

c 

■*- 

« 

• 

#— 

ft 

ft 

X 

*** 

2 

< 

< 

< 

X 

a 

►— 

tt 

X 

w 

IT 

tt 

t- 

-J  z 

LLj  l. 

C  " 

x 

t- 

Q  * 

<  II 

X  C 
C  >. 

if 


U  I*  ►-  n 
C  -  a  - 
«  .  c  X 

—  p-  if  — 
ir  i 

•  *  •  Q, 

I  •  I' 

X  C  '<  n 

> <  >w  p*t  p*> 


<  l~ 

—  a 


a  ii  z 
v 


c 

* 


c 

X 

« 


<v 


if  e 

«  —  *-  c  c 


a 

z 

< 

« 

X 

t- 

u 

« 

z 

If 

« 


•—  I 
X  » 

« 


u. 

c 

# 


It  X  * 
ff  —  - 
<  —  <\ 
i_  a  a 


w 

C 

cc 

X 

z 

a.  Z 

UN 

• 

X 

«  a 

rnm 

z 

IT 

*■> 

tt 

i/* 

w 

I  * 

• 

a 

• 

a 

» 

• 

CL 

c 

n 

11 

a 

* 

H*  pM 

c 

1 

c 

c  ^ 

X 

p- 

C 

c  a 

““ 

z 

IT 

*— 

II 

1 

> 

tt 

<  c 

u. 

It 

1 

z 

c 

•I 

X 

X 

2 

H* 

a  *- 

•— 

• 

#-> 

tr  !>• 

c 

X 

X 

1 

1 

*"*■ 

w 

•— 

1- 

< 

* 

X 

c 

H  II 

ii 

X 

x 

X 

X 

C 

u 

2 

X 

• 

•“  c: 

!»* 

X 

w 

X 

X 

< 

1/ 

Z 

i— - 

tt 

« 

<1 

tt 

>— 

-  Z 

•_ 

tt 

1! 

It 

II 

*• 

2 

< 

c 

a 

c 

_i  c 

mm 

#— 

c 

mm 

^  •— 

c 

X 

X 

<\ 

a 

< 

tt 

c 

0 

a 

-  c 

X 

X 

c 

<c 

fv  */ 

a 

c 

X 

X 

cc 

a 

l_ 

» 

IN 

X 


x  <■ 

—  H  _■ 
Ui  J  U.1 

>  U.  c 

«  C  I 
u  ♦ 

V  X 
-J  X  X 
•I  II  II 
<  *—  r\ 
c  x  a 


z 

-  c 

cc  i— 

»  h- 

(\  — 

*  -J 

-  C  - 

*  if  y 

t“  cv 

Z  H  <— 
•—  u.  • 

*  C  a 

X  » 

=  C  y 

C  K  • 
<  < 

»  X  • 

is  Z 

a  a  . 

'  •-  2 
*-  <  _ 
cr  zo 

U- 

U  U.  2 

>  u-  X  •— 
<  " 

^  2  C 

X- 

x-  Of  I 

-  2  a  _ 
<  C  >  <  ! 
c  c  z  e  i 


<* 


c 


81 


1363-6-T 


u  « 
c  "» 
►  — 
z  < 

1/1  ft 

°S 

•  c 

i/l  x. 
1/  X 
.  Z 
i/.  r 
o  u. 
►  » 
_  >. 
- - ’  w 

a  o 

•  V 
—  m 

a  < 


st  < 
x  C 
>  s 
—  a- 
x  l 


< 

— 

— 

«k 

>K 

• 

X 

X 

X 

CsJ 

y 

u> 

Ik 

u. 

A 

< 

ft 

c 

c 

C 

e 

h- 

ft- 

z 

z 

Z 

z 

< 

* 

— • 

w 

»»— 

- 

C 

< 

• 

» 

et 

U. 

«Si 

* 

Z 

z 

3 

ft 

ft 

ft 

• 

< 

ft 

V 

»- 

h- 

*~ 

ft 

K 

c 

• 

» 

• 

z 

Uu 

if 

< 

2 

z 

z 

• 

1 

O 

** 

z 

H 

»■ 

• 

a 

a 

a 

X* 

• 

c 

rsj 

• 

♦ 

• 

r 

•x 

41 

h 

►* 

H 

c 

c 

2 

a 

z 

2 

2 

L 

fv 

• 

• 

# 

to 

to 

to¬ 

• 

«** 

a 

c 

to 

• 

• 

ft 

Cl 

X 

2 

z 

< 

> 

> 

> 

— -  N 

IN.  I** 

* 

a 

• 

<s 

K 

■* 

<«  * 

ft 

• 

to 

«► 

1/ 

• 

c 

c 

c  c 

►-  * 

ft 

< 

C 

z 

9 

X 

4| 

48 

< 

►-  * 

♦  ♦ 

♦ 

< 

♦ 

V 

Z 

• 

• 

ft 

0 + 

**  to 

IV  *> 

— 

IV 

* 

c 

«/ 

ft  C  C 

to»  to 

«  - 

— » 

« 

u. 

c 

X 

v* 

4t 

cr 

* 

a 

P- 

c  c 

4  t 

ft  ft 

1 

K 

P 

A 

>* 

• 

• 

• 

ft 

► 

♦ 

*  * 

♦  ♦ 

2 

4 

< 

w 

X 

> 

to 

to 

to 

to 

«—  to 

•  * 

* 

< 

< 

K 

c  z 

or 

tf 

« 

y 

-■  z 

to*  ft— 

«v  * 

*» 

4/ 

C  X 

■** 

*■" 

• 

to*  to* 

*  • 

ti- 

a 

z 

< 

• 

•  it 

%r 

i/*. 

•  • 

i—  — 

ft 

ft1- 

Z 

c 

c 

X 

< 

a 

a 

Cl 

to* 

o  c 

Cl  9 

« 

to* 

to 

'V 

ft- 

X 

z 

X 

4 

U.  U- 

to  *to 

«  ft 

— 

to* 

O  '« 

4C 

c 

ax 

to* 

ax 

tt 

• 

ft 

to- 

*  if 

a 

N  #* 

•to 

a 

E 

z 

z 

•  u. 

z 

c  z 

40 

t1 

K 

t- 

c 

u 

to 

♦ 

i 

to 

to*- 

» 

to 

z 

It 

c 

</ 

a 

ft 

« 

a 

z 

Z 

z 

X 

X 

X 

X 

x  C 

z  c  z 

to 

« 

« 

X 

X 

✓  0 

«L 

ft 

* 

> 

P 

p 

-J 

a 

t! 

— J 

ft 

to- 

to* 

ft 

*  ft 

ft 

»  ft 

ft 

• 

X 

«r 

to- 

X 

c  x 

• 

P 

ft 

< 

_j  ft 

41 

-J 

ft- 

ft- 

ft* 

•X 

ft* 

ft* 

ft- 

p* 

c 

c 

c 

•  3  C 

to- 

fv. 

c 

z 

a 

« 

u 

49 

it 

ft 

to* 

u. 

u. 

to* 

-  C 

ft* 

-  c 

ft* 

ft* 

u 

Z 

i/ 

c 

•  t 

u 

a 

a 

< 

z 

c 

z 

C 

z 

e 

< 

to 

** 

< 

ft  c 

m 

ft  c 

< 

*C 

a 

u. 

—  fv 


83 


j 


1363-6-T 


* 

z 

QC 

< 

f*  # 

z 

KJ  2 

X. 

!*>  C 

1/) 

#  »/! 

c 

z  + 

, 

c  — 

#• 

X 

z 

•f  (- 

a 

—  ur 

< 

X  * 

*— 

K  X 

ur 

X 

C 

♦  1 

u 

X  _ 

o 

X  * 

4 

—  z 

** 

■»  LL 

z 

z  c 

c 

u.  * 

X 

X 

ur 

c  - 

X 

X 

4 

«  • 

* 

« 

Z 

^  r" 

c 

c 

LL 

*  c 

•  W 

C  - 
— 

—  rv 
■tt  rv  »~ 


X  X  c 
■v  v  I 

a  a  * 
*  *  - 
—  —  a 
,  •  a 

*-  f  < 


rv 

9 

* 

*** 

fv 

V 

X 

• 

• 

z 

• 

fv 

— 

X 

c 

c 

X— 

<■** 

V. 

CV 

««w» 

w 

■*— ■ 

i/: 

1 

c 

— * 

fvl 

a 

4 

4 

C* 

♦ 

r— 

> 

X 

«— 

'v 

• 

► 

fv 

LL 

• 

• 

— ■ 

c 

• 

> 

a 

C 

c 

C 

a 

c 

w 

X 

tt 

» 

» 

a 

«*• 

— 

a 

u 

• 

• 

< 

a 

tt 

X 

«— ■ 

t*4 

*- 

»— • 

X 

*-* 

a. 

z 

«•» 

w 

i/* 

a 

c 

c 

4 

f  f 

tt 

tt 

C 

c' 

\J 

-«• 

— 

1/ 

*-* 

c 

tt 

a 

A 

+ 

*— 

> 

c 

i 

*- 

z 

2 

X 

1 

X 

r- 

LL 

a 

i 

*- 

h- 

sS 

X 

c 

c 

> 

X 

U. 

U' 

V— 

<\ 

f~- 

X 

1 

>>+ 

w 

IT 

4 

X 

N4  fv' 

1 

z 

tt 

tt 

tt 

• 

fC* 

“ 

X 

fv! 

rv 

> 

<~ 

c 

a  »r 

c 

c 

c 

,»  LL 

h- 

X 

X 

tt 

# 

> 

a 

</ 

*•  X 

a- 

X 

X 

» 

^  ♦ 

u 

« 

■tt 

CV 

(N 

X 

w 

i  a 

"v 

V 

• 

2  ^ 

e 

X 

X 

> 

X 

a 

il 

tt 

—  4 

2 

*•— 

*— 

c 

C  u- 

u 

X 

X 

> 

X 

X 

c 

LL 

rr 

*— 

a 

a 

w 

t/“  w 

II 

II 

II 

11 

II 

u 

tt 

h* 

< ^  LL 

X 

n 

it 

II 

4  •• 

f\ 

•“ 

f“ 

r- 

fv 

1! 

r<" 

1 

fv  H 

z 

a 

2 

fT 

2  5 

K 

X 

X 

fv 

fv 

LL 

a 

*V  fT 

c 

a 

a 

K 

LL  2 

K 

X 

X 

> 

> 

H 

a 

LL 

c 

< 

< 

a 

C  u 

l\ 

U. 

*-  c 

a  2 


86 


<;impni)TINF  C  T  MPC  (  s,b,f,MAv]  .  YINTH.XIP,  I IMDFK  ) 
Twt<;  <;tmRniiT!NP  if  fop  t no  f t;p m ] hn 

r.  H  W  P  L  F  X  FAFtA*XT  P  ♦  X .!  ,  FNFWX,  Xf  •  F  ,  X  I  M  T  G 
P  APR=X  1  MTT,  (  A  )  4-x  I  MTT,  <  P,  ) 


1363-6-T 


rr 

r— 

— 

c 

<£ 

tr 

c 

c 

C 

cc 

t- 

C 

< 

F“ 

c 

c 

c 

►* 

c 

c 

• 

K 

tx 

c 

• 

•m 

mJ 

e 

C 

* 

# 

mm 

c 

• 

#«• 

X 

X 

rr  — 

3 

IT 

< 

— *  m 

3  C 

a 

X 

—  X 

a 

2 

• 

a 

tr  1 

2  • 

IX 

* 

• 

*  a 

X  K 

4 

x  H 

• 

•  »_ 

-  C 

tf 

• 

4  C 

tr  x 

«v  X 

a:  c 

u.  a  •  *- 

tr  «  ►-  rr  2 

•  u.  x  v  •— 

•»  «  4  tx  C  X 

<  l  a  <  -  ii  ii 

|  X  I-  II  X  X  X 

I  II  x  .111  a  J 

n  a  II  U.  «  C  IX 

X  I - .2X2  2 

X  >  X  X  X  *-  u. 


|f  X  I  •  mi  + 

o  a  «  •  X 

»-  c  x  x  —  x  tr 

x  z  u;  it  x  •  3  4  ££ 

+  *-  C  »  4  C  U.  ~  < 
— .  II  2  *  <  *  Z  X  tj 

X  X  •-  I  II  •  x  ~ 

«.  u  «  x  3  c  ti¬ 
ll  C  n  a  —  -  cl 
-,2a  X  2  J'u-i-u. 
x  •—  >-  x  x  tr.  —  x  — 


X  . 

x 

*  +  tr 

C  X  x  c 

C  3  IXC 

^  u-  r  *- 
&  z  x 

r-  C  tx  it  c 

X  h  4  J  I" 

ii  tr  a 

-  C  H  2  C 

x  e  tr  x  c 


>c  f\. 

c  c  c 
c  c  c 


if  h  ec 

c  c  c 

c  c  c 


88 


f.FT  <;t!R=r.niUMN  OF  7FPn<;  IN  THF  K=TH  COLUMN 
70  if  f  K  .FO.  N>  r,n  TO  1 


Z  Z  -  - 


*  X 
It  It  It  — 

—  a. 

1  II 
r-  or  »  — 

c  c  -  ~ 
C  C  *  cr 


z 

c 

*■« 

K 

Ml 

*- 

• 

X 

K 

i 

*** 

IT 

1 

< 

cr 

** 

¥ 

cr 

«Ci 

rv 

» 

cc 

* 

# 

-5 

*— 

* 

* 

• 

*— 

<VJ 

< 

* 

a 

* 

W 

1 

♦ 

#— 

< 

M— 

« 

> 

2  * 

X 

w 

c 

cc 

• 

» 

< 

» 

*■* 

C  "5 

cr 

»— 

1 

• 

*  ar 

•  za-. 

i 

w 

c 

ii  » 

C  1 

X  X 

►— 

Ui  x  r 

» 

II  ♦ 

Uj  ^  Ul 

•I  — 

-  2  * 

X  —  I 
•  H  X 

—  7 

-  C  7  C 
«  u  —  c 

« 

r-  h 
a 

c 


H  •+  f 

*“  X 
C  •  II 
7  C 


ii  a  -  a 
x  x  v  — 


\f  w  x  —  *“ 

I,  II  •-  II  —  — 

lCTII'»«X*"*'HIe  1 
Ca  2  Z  it  2  m  of  I'  T 

—  x.  O  \r  it  II  «v&  —  oc 

II  _  -  II  F  ^  I*  I  1 - - - 

Z  X  7  X  7  X  -_H 


0 c 
< 
o 

z 

c 


c 

h 

c 

c 


c 

• 

o 

a1 


c 

41 

z 


c 

7 

< 

c 

* 

c 

a  ~ 
•  xi 


co 


c 

c 

CM 


< 

I*- 

CV 

IT 


c 

c 

4 

rv 


=  c  -  c  c 
v  c  x  c  a  c.  cr  c 
a 


►■Pc 

».  C  •-  -  a  Z 
x  a  u. 


c  < 

*-  IV 


c 

cr 

c 

CJ 


X  X:  •/. 

«J  Xj  —  tC 

Z  _!  <  C  • 

2  X  <  O  «V  c 

u.  u.  U.  II  HZ 

jut  i-  Cx  b 
_ia~i/l-tr- 

<  X  *  tf  Pc 

aCa«C<aZ 
o  t  «-  t  c  c  a  ul 

c  c 


c 

c 


< 

c 

* 


doc  m 

IS  NT  CONTROL  DATA  ■  R 

4D 

The  University  of  Michigan  Radiation  Laboratory,  Dept,  of 
Electrical  Engineering,  201  Catherine  Street, 

Ann  Arbor,  Michigan  48108  _ _ 


3.  REPORT  TITLE 

SCATTERING  BY  A  PERIODIC  SURFACE 


UNCLASSIFIED 


4.  OBtCRlRTivc  noth  (Type  ol  t*poH  md  Inclusive 

Scientific  Interim  _ 


e.  AUTHORlSl  (ritii  tunnb,  middle  Initial,  laal  name) 

Tommy  C.  Tong 
Thomas  B.  A.  Senior 


«.  report  oat* 

September  1970 


8«.  CONTRACT  ON  GRANT  NO. 

F19268-6 8-00071 

b.  PROJlC  T  NO.  ,  Task  No. and  Work  Unit  No. 
5635-02-01 

c.  DoD  Element:  61102F 

„  DoD  Subelement  681305 


»«.  TOTAL  NO.  OF  PACE* 

96 


lb.  NO.  OF  ftCFt 

30 


M.  ORISINATOn'S  REPORT  NUMBERISI 

1363-6-T 

Scientific  Report  No.  7 

<th.  OTHER  REPORT  NOisi  (Any  other  number*  thot  nuiy  be  maalfrutd 
f/i/a  report) 

AFCRL-70-0507 


10.  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution  is  unlimited. 


II.  SUPPLEMENTARY  NOTES 

TECH,  OTHER 


12  SPONSORING  MILITARY  ACTIVITY 

Air  Force  Cambridge  Research  Laboratories(CR 
Laurence  G.  Hanscom  Field 
Bedford,  Massachusetts,  01730 


It.  ABSTRACT 


Numerical  procedures  are  developed  for  the  digital  solution  of  the  integral 
equations  for  the  currents  induced  in  a  perfectly  conducting,  two  dimensional 
periodic  surface  when  a  plane  electromagnetic  wave  is  incident.  Data  are 
obtained  for  both  the  surface  and  far  fields  for  a  variety  of  sinusoidal  surfaces 
for  plane  waves  of  either  polarization  at  oblique  as  well  as  normal  incidence, 
and  the  results  are  compared  with  the  predictions  of  physical  optics. 


UNCLASSIFIED 

curUy  Cltiimtflollon 


1363-6-T 


Periodic  Surfaces 
Diffraction  Gratings 
Surface  Currents 
Diffraction  Modes 
Numerical  Techniques 
Exact  Data 
Physical  Optics 


